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Abstract 



In this paper we study the Yang-Baxter integrable structure of Conformal Field Theories 
with extended conformal symmetry generated by the W3 algebra. We explicitly construct 
various T and Q-operators which act in the irreducible highest weight modules of the 
W3 algebra. These operators can be viewed as continuous field theory analogues of the 
commuting transfer matrices and Q-matrices of the integrable lattice systems associated 
with the quantum algebra U q (sl(3)). We formulate several conjectures detailing certain 
analytic characteristics of the Q-operators and propose exact asymptotic expansions of the 
T and Q-operators at large values of the spectral parameter. We show, in particular, that 
the asymptotic expansion of the T-operators generates an infinite set of local integrals of 
motion of the W 3 CFT which in the classical limit reproduces an infinite set of conserved 
Hamiltonians associated with the classical Boussinesq equation. We further study the 
vacuum eigenvalues of the Q-operators (corresponding to the highest weight vector of 
the W3 module) and show that they are simply related to the expectation values of the 
boundary exponential fields in the non-equilibrium boundary affine Toda field theory with 
zero bulk mass. 
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1 Introduction 



The Yang-Baxter equation and its associated integrability structures play an important 
role in the understanding of Integrable Quantum Field Theory (IQFT) |T]. Recently, the 
authors of [0-0] developed a continuous version of the Quantum Inverse Problem Method 
(QIPM) H or Baxter's commuting transfer matrix method [^] to study IQFT, especially 
the finite volume IQFT. This new method inherits all the analytic power of the theory 
of integrable lattice models but does not require intermediate lattice regularization and 
allows from the very beginning to work in terms of the continuous IQFT. The method has 
been most successfully applied [|2|-|J to the c < 1 Conformal Field Theory (CFT) which 
is a simplest nontrivial example of IQFT. The fundamental problem addressed in these 
papers is the simultaneous diagonalization of an infinite set of mutually commuting local 
integrals of motion (IM) in CFT. 

In this paper we extend the results of p|-§J to CFT's with extended conformal sym- 
metry generated by the simplest P^-algebra W3 ||. Such a CFT involves additional local 
spin-3 fields W{u) and W(u) besides the usual left and right components T{u) and T(u) of 
the stress-energy tensor. In what follows we will concentrate on the left chiral component 
W3 of the full symmetry algebra W3 x W3. The algebra W3 is generated by the fields 
T{u) and W (u) along with all composite fields built from products of T(u) and W(u) 
and their derivatives. We will interpret the variable u as a complex coordinate on a 2D 
cylinder of circumference 2n and impose the periodic boundary conditions 

T(u + 2tt) =T(u), W(u + 2ti) = W{u). (1.1) 



The algebra W3 possesses an infinite-dimensional Abelian subalgebra ||,[T(J generated 
by the local integrals of motion (IM), I*. G W3 

Ifc = / — T k+1 (u), keZ + , k^O (mod 3) (1.2) 
Jo 27r 

where the densities T^{u) are appropriately regularized polynomials in T(u) and W(u) 
and their derivatives. The first two densities read 

T 2 (u) = T(u), T 3 (u) = W(u), (1.3) 

while the higher densities are determined by the requirement of commutativity 

[Mi] = 0, (1.4) 
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and by the "spin" assignment |§ (see Section || for the details) . 

We explicitly construct operator valued functions T(t) and T(t) where t is a complex 
parameter. These operators act invariantly in each irreducible highest weight module of 
the W3 algebra 

T(t),T(t): V A2iA3 - V A2iA3 (1.5) 

where A 2 , A 3 are the highest weights and they commute among themselves for any values 
of t, i.e., 

[T(t),T(t')] = [T(t),T(t')} = [T(t),T(t')] = 0. (1.6) 

These operators are defined in terms of certain monodromy matrices associated with the 
three-dimensional representation of the quantum algebra U q (sl(3)) where 

q = e*", (1.7) 

and g is related to the central charge c of the W3 algebra as 

c = 50-2A(g + g- 1 ), (1.8) 

Here we mostly restrict our considerations to the "quasiclassical domain" 

-oo<c<-2, (1.9) 

however, the results are valid in a wider region — 00 < c < 2, corresponding to < g < 1. 

Obviously, the operators T(t) and T(t) are W3 analogues of Baxter's commuting trans- 
fer matrices from lattice theory. They enjoy remarkable analyticity properties, namely, 
they are entire functions of t with an essential singularity at t = 00 and their asymptotic 
expansion near this point is described in terms of the local IM ( |1 . 2| ) . For instance, 

00 

logT(t) ~mt^) +X! C fc t~^j (1.10) 

fc=i 

where m and C' k are some explicitly known coefficients which depend on the central charge 
c only and C' k = for k = (mod 3). A similar expansion holds for T(t). Therefore the 
operators T(t) and T(t) can be thought of as generating functions of the local IM since 
their eigenvalues contain all the information about the eigenvalues of the local IM. 

Another powerful method known in the lattice theory is based on the so-called Q- 
operators. This method was introduced by Baxter in his original solution of the eight- 
vertex model of statistical mechanics [|11|]. Here we construct appropriate W3 versions of 
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Baxter's Q-operators (actually, there are six such operators Qi(t) and Q»(i), i = 1,2,3). 
They are defined again as traces of certain monodromy matrices this time associated with 
infinite-dimensional representations of the so-called g-oscillator algebra. The Q-operators 
satisfy the following relations with the operators T(t) and T(£) discussed above 



Q l (tq 3 )-T(tq 1 2)Q i (tq) + T(tq--2)Q i (t q - 1 )-q i (tq- 3 ) = 0, 



Q i (V)-T(*g5)Q i (tg)+T(tg-5)Q,(t g -i)-Q i (ig-^) = 



;i.n) 



where % = 1,2, 3. These relations are obviously the {7 9 (sZ(3))-analogues |12|,|T3| of Baxter's 



famous T-Q relation. The relations are just the first representatives of a rich 

family of the functional relations obeyed by the Q-operators in the U q (sl(3)) related 
integrable systems. In particular, in Section 4 we define more general T-operators T M (i) 
and T M (£) corresponding to finite-dimensional irreducible representations of U q (gl(3)). 
Here \i = /i 2 , Ate) is the highest weight of the representation such that /xi — /x 2 and 
/i 2 ~ Ate are non- negative integers^. All these T-operators can be expressed through Qj(t) 
and Qj(t) by the means of only two elegant relations 



T (W)/12j/13) (t) = z^det || Qi(tg 



V^2,w)W = V det II Qi(*?' 



+2/Lt' 



-2/i' 



Uj=l,2,3 

(1.12) 

Lj=l,2,3 



where fj,\ — /ii + 1, /i' 2 = /i 2 , /i' 3 = /i 3 — 1, and where z is a simple operator acting 
diagonally in each W3 module which can be interpreted as the "quantum Wronskian" |3| 
of the third order difference equations (|1.11|) 



z = det 1^(^)11,^3 = -det \\Qi(tq 2j )\\ i>j=1)2j3 (1.13) 



We prove the relations ( 1.12j) algebraically by using decomposition properties of direct 



products of certain "g-oscillator" representations of the Borel subalgebra of U g (sl(3))f\. 
These relations provide a very simply tool for derivation of various "fusion relations" 

1 The existence of the two different sets of operators T M (t) and T^(t) is related to the fact that there 
are two non-equivalent "evaluation maps" from U q (sl(3)) to U q (gl(3 j). 

2 It should noted that some particular cases of (|l,12|) were previously obtained Jl^,|l^] from the prop- 
erties of the corresponding eigenvalues determined by the so-called "nested Bethe Ansatz" solutions [151 . 
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for the T-operators |16|,|l7|. In particular, the so-called RSOS-reductions |L6| of the T- 
operators simply follow from ( 1.12| ) when two columns of the three-by-three matrices 
therein become proportional. 

The eigenvalues of the Q-operators can be studied by various methods. In particular, 
they satisfy the Destri-de-Vega non-linear integral equations fllSl-pOfl. We solve these 
equations exactly in a particular case when A 2 , A 3 — > oo, where A 2 , A 3 are the highest 
weight in ( |1.5|) . On the basis of this result in Section 6 we conjecture exact asymptotic 
asymptotic expansions of the eigenvalues of the Q-operators at large values of t. These 
expansions exhibit remarkable duality properties of the Q-operators with respect to the 
substitution g — > g^ 1 similar to the case of c < 1 CFT ||. We support our conjectures 
by consideration of certain exactly solvable cases and numerical calculation of the zeroes 
of eigenvalues of the Q-operators. 

As shown in |2| the T and Q-operators define integrable boundary conditions [^1] for 
CFT in the geometry of the half-infinite cylinder. In particular, the vacuum eigenvalues of 
these operators can be interpreted || as finite-temperature partition functions of certain 
integrable QFT models with a non-trivial boundary interaction. These models are of 
interest in many branches of physics. For example, the integrable boundary sine-Gordon 
model pil] (with zero bulk mass) finds important application in dissipative quantum 



mechanics p^p5|; as explained in [g6],|27] it also describes the universal current through 
the point contact in the quantum Hall system [p|,^8|-0O|. Here we consider a simplest 
generalisation of this model, namely, the boundary affine Toda theory with zero bulk 
mass. Its Hamiltonian has the form 

H= J- f rfx(n 2 + (^) 2 )-f £ yj e^>*W, (1.14) 
^9 J -oo ^9 =1 

where g, k are parameters and $(x) = ($i(x), $ 2 (x)), n(x) = (Ili(x), U 2 (x)) are field 
operators obeying canonical commutation relations 

[n a (x),$ 6 (x')] = -2mg 5{x - x')S ab , a, b = 1,2. (1.15) 

and Cj, j = 1,2,3 are the simple roots of the algebra A% ■ The "amplitudes" yi,y 2} y3 in 
(|l-14j) are not c-numbers, but quantum operators describing some additional boundary 
degrees of freedom. They satisfy the commutation relations 

3/13/2 = 93/23/1, 2/22/3 = <72/32/2, 2/32/1 = 93/12/3, 3/13/23/3 = 1, (1.16) 

where q = e mg . It should be noted that the integrability of the boundary affine Toda sys- 
tems has been studied previously both on classical |J]J and quantum levels [32]. In |JT 
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it was shown that (with the exception of the boundary sine-Gordon model) the integra- 
bility in the boundary affine Toda theory imposes severe restrictions on the amplitudes 
of the boundary interaction. Namely, the boundary interaction should vanish when the 
bulk mass is zero. This conclusion was then confirmed |32) with the perturbation theory 
calculations in the quantum case. We would like to stress here that the model ( |1.14|) does 
not fall under this "no-go theorem" since ( |1.14| ) involves additional boundary interaction, 
which was not considered in [31], j32|. We claim here that the theory (|1.14| ) is integrable, 
but postpone detailed considerations of this question to a separate publication. Actually, 
the sufficient condition for the integrability of ( |1.14| ) is that operators t/j obey the Serre 
relations of quantum affine algebra U q (A^) 

Vi Vj ~(q + l" 1 ) Vi Vj Vi + Vj Vi =0, i # j. (1.17) 

and not necessarily restricted to the case (|1.16|) . A Similar statement holds for general- 
isations of ( |1.14Tj to other affine algebras. In Section |7.1| we consider a non-equilibrium 
variant of Ql.14 ) with a particular time-dependent boundary interaction. We show that 
expectation values of boundary exponential fields over a stationary non-equilibrium state 
in this theory are simply related to the vacuum eigenvalues of the Q-operators of the 
W3 CFT. This calculation generalises the corresponding result of for the boundary 
sine-Gordon model. 



2 Quantum Boussinesq theory 

2.1 Extended W3 conformal symmetry and local integrals of mo- 
tion 

In this paper we consider CFT with extended conformal symmetry generated by the 



simplest W- algebra W3 (sometimes called WA 2 ) f|,|33]]. Such a CFT involves additional 
local spin-3 fields W{u) and W(u) besides the usual left and right components T{u) and 
T(u) of the stress-energy tensor. The left chiral component W3 of the full symmetry 
algebra W3 x W3 is generated by the fields T{u) and W{u) along with all composite fields 
built from products of T(u) and W(u) and their derivatives. We will interpret the variable 
u as a complex coordinate on a 2D cylinder of circumference 2ir and impose the periodic 
boundary conditions ( |1.1|) . The Fourier mode expansions of T{u) and W(u) 

+00 +00 
T(u) = + E L -n emZ > W ( u ) = E W - emZ > 
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are expressed in terms of the operators L n and W n satisfying the commutation relations 



[L n , L m ] 
[L n ,W m ] = (2n - m) W n+ 

[W n ,W m ] 



(n - m) L n+m + — (n 3 - n) 5 n+mfi , 



1 , (n — m) 

-(n - m)A n+TO + — ^ — 
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--(n + 2)(m + 2) 
6 



+ 



[n + m + 2) (n + m + 3) 



(n 2 - 4)(n 2 - l)nS n+m ,o, 



(2.2) 
(2.3) 



1080 6 2 



where c is the central charge, 



b 2 = 16/(22 + 5c), 



and the operators A n are expressed through L n as 

+00 ^ 

A n = ^ ] : : +— x n L n , 

with 



(2.4) 



fc=— oo 



z 2 , = (1 + /)(1 - I), x 2l+1 = (2 + 0(1 - 0- 



The symbol : : in ( |2.4|) denotes the normal ordering such that operators L n with bigger n 
are placed to the right. 

The chiral space of state TC c hirai is built up from a suitable collection of highest weight 
modules 

7~Lchiral = ® H A («) . (a),, (2.5) 

where the parameters A 2 and A 3 are the highest weights, and the associated highest 
weight vectors |A 2 , A 3 ) G V(a 2 ,a 3 ) satisfy the equations 

L n \ A 2 , A 3 ) = 0, W n \ A 2 , A 3 ) = 0, n > 0; (2.6) 

L |A 2 ,A 3 ) = A 2 |A 2 ,A 3 ), iy |A 2 ,A 3 ) = A 3 |A 2 ,A 3 ). (2.7) 

The algebra W 3 possesses an infinite Abelian subalgebra H|l(| generated by the local 
integrals of motion (IM), G W 3 

I* = 



— T k+1 (u), ke z_, 

Z7T 



fc ^ (mod 3) 



(2.8) 



where the densities Tk(u) are appropriately regularized polynomials in T(u) and W(u) 
and their derivatives. The integer k + 1 in ( |2.8| ) denotes the spin of the local field T k +i 

dw 



2-ni 



(w - u)T(T(w)T k+1 (u)) = (k + 1) r fc+1 (u), 



(2.9) 



where T stands for the chronological ordering 



n-t u lP , u . A(w)B(u) , if %m u > %m w ; 
T(A{w)B{u)) = { (2.10) 
B[u)A[w) , if Xsm w > <sm u . 



The set of possible values of k in ( |2.8| ) is determined by the allowed values of spins of the 
densities Tk+i(u) 

spin[T fc+1 (M)] = k + 1, k G Z + , fc ^ 3 (mod 3) (2.11) 

The first few densities T^iu) can be written as [ID|^| 



T 2 (u)=T(u), T 3 (u) = W(u), T 5 (u)=:T(u)W{u):, (2.12) 

r — 1 n 

T 6 (u) =: T( M ) 3 : +9 : W(u) 2 : +-^- : T\uf :, (2.13) 

where the prime denotes the derivative and : : denotes the regularized product of the fields, 
such that the product : A{u)B{u) : is defined as the regular term (of order 0((u — v )°) in 
the operator product expansion of A(w)S(t)).Q The densities ( |2.13|) are defined up to total 
derivatives which, obviously, have no effect in fl2.8|) . A general expression for the densities 
( p,13| ) is not known, however, once the first two densities in ( |2.13| ) are fixed all other 
higher densities can be uniquely determined (up to a normalisation) by the requirement 
of commutativity 

[Ik, M =0, (2.14) 

and the spin assignment fl2.11|) . Alternatively one can formulate the last requirement as 
the "degree assignment", namely, if one assigns degrees 2 and 3 to the fields T(u) and 
W(u), respectively, and the degree 1 to each derivative then T fc+1 (w) is a homogeneous 
polynomial of the degree k + 1. We normalise the densities Tk+i(u) by fixing coefficients 
to the terms containing highest power of T(u), 

T 2n (u) = (T(u)) n + ..., T 2n+1 (u) = (T(u)) n - 1 W(u) + ... J n = 1,2,3,... (2.15) 



The integrals ( p.8|) define operators Ik- Va 2 ,a 3 —* Va 2 ,a 3 which can be expressed in 
terms of the generators L n and W n , for example 



3 A misprint in the expression for Tq(u) in [T^] is corrected here. 

4 We assume that these regularized product are calculated in the cylinder geometry, as they are not 
invariant under conformal transformations. 
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h 
I 4 



u 

Wo 

+00 



24 



J2(L~nW n + W. n L n ) + W L - °-^W , 



n=l 



24 



;2.16) 



/ 14 + c 2 c\ -3 

■ L ni L n2 L n3 : + 2 I — — n — 4 J 2 ' -^l-2r-^2r 

™i+™2+™3=0 n=l ^ ^ r=l 

+l8 f „ + 9W? _ £+8^ + (£±i^l io _ e(e+23)(7e + 30) 
L^i 8 1Q9 96768 



71=1 



Here the symbol : : denotes the same normal ordering as in ( |2.4j ) (where L n with bigger 
n are placed to the right). Note that although these operators are not polynomials in L n 
and W n their action in Va 2 ,a 3 is well defined. 



2.2 Relation to classical Boussinesq theory 

In this paper we will study the problem of simultaneous diagonalization of the infinite set 
of the operators I*, in V(a 2 ,a 3 )- This problem can be regarded as the quantum version of 
the related theory of the classical Boussinesq equation ]34j , which is a well known example 
of a completely integrable Hamiltonian system possessing an infinite set of local IM [|35] . 
As is known || the W-algebras can regarded as the quantization of the Gelfand-Dickey 
Poisson bracket algebras |36|]. The classical limit of (|2.2j ) corresponds to c — ► — 00. In this 
limit the substitution 



T(u) 



24 



U(u), W {u)^i(-^fv{u), 



48vr 



ic 



where 



U{u + 2tt) = U{u), V(u + 2n) = V{u), 



(2.17) 



(2.18) 



reduces the algebra (|2.2|) to the following Poisson bracket algebra j37 



{U(u),U(v)} 
{U(u),V(v)} 
{V(u),V(v)} 



+ 



-(U(u) + U(v)) 5'(u -v) + 2 5"'(u - v), 
-(V(u) + 2V(v)) 5'{u-v), 
\{U"{u) + U"{v)) + l -{U\u) + U\v))] 8'{u 

^-{U{u) + U{v))5"'{u -v)- U^\u - v), 
12 6 



(2.19) 
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which is a particular example of the Gelfand-Dickey bracket algebra. The substitution 
( 2.17|) brings the local IM ( |2.8|) to their classical counterparts 



I^r^f-^^/f , c--oc (2.20) 



which form a commutative Poisson bracket algebra {/^ , ij } = . For example, the 



first few classical IM 7^ read 



T (cl) 



1 2 



\ ^-U(u), (2.21) 
Jo ^ 

rp-V(u), (2.22) 
Jo ^ 

4 C ° = r^U(u)V(u), (2.23) 

U{uf + W{uf + V(u) 2 ). (2.24) 



T {d) 



The Poisson bracket algebra ( |2.19| ) is known to describe the second Hamiltonian structure 
of the Boussinesq equation [[Tj],|38| provided one of the infinite set of the classical IM Ijf' 
is taken as a Hamiltonian. 
The evolution equations^ 

\U = {li d \U}, d tk V = {l£ l) ,V}, t 1 = u, k^O (mod 3), (2.25) 

describe isospectral deformations of the third-order differential operator 

L = dl - U(u) d u + V(u) - Uj'(u) - A 3 , (2.26) 

where A is the spectral parameter. In particular, if one defines the monodromy matrix 
M(A), which is an element of the group SL(3), as 

(V>i(u+27r),<0 2 (u + 27r),^(« + 27r)) = (^i(w), ^(u))M(A), (2.27) 

where ■0i(w),^ 2 ( w ) an d ^(u) are three linearly independent solutions of the equation 
Lip = 0, then the eigenvalues of M(A) are in involution among themselves and with all 



The original Boussinesq equation of 



d t U = -2d u V, d t V = ^d u 



\d 2 u U-U 2 
2 " 



corresponds to k = 2 in (2.25) 
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local IM ( p.24|) (with respect to the Poisson bracket algebra ( |2.19|) ). The trace of the 
monodromy matrix 

T (d) (A) = Tr M(A), (2.28) 

can be regarded as a generating function for the local IM (|2.24| ) as it expands in the 
asymptotic series |38| 



f logT^(A) = A + ^c fc 4 c/) A^ 

k=l 



A — > +00 



;2.29) 



where (c^ = for k = (mod 3)) are some numerical coefficients. The simplest 
way to calculate these coefficients is to consider the case when both U(u) and V(u) are 
constant (i.e., independent of u) and take real values. In this case ( |2.29| ) reduces to the 
A — > +00 expansion of the positive root of the cubic equation x 3 + 2Ux + V — A 3 = 
which can be easily calculated. Similarly to ( p,15|) we normalise the classical IM 7^ by 
the term containing the highest power of U(u) 



r(cl) 



2~ 



du 



< 2 „-. ■ / 7r(u(u)r+. 

/ 27T 



j{d) 



2n 



du 
2^ 



(U(u)T- 1 V(u) + ... 



n = 1,2, 3, 



(2.30) 

Note that this normalisation is consistent with (|2.15|) , ( p. 171 ) and ( |2.20|) . It follows then 
that 



C-2n-l 



4vrcos ((i^h)r(2^i; 



3n\T 



2-n-- 



C2n+2 



4?r cos(ggg)rgg^ 

n!T(-^±i) 



(2.31) 



where n = 0,1,2,.... The required calculations are given in Appendix C; the above 
expressions are simple corollaries of ( |C.9| ). 

For our purposes it will be more convenient to work with a first order differential op- 
erator instead of fl2.26|) . To define it we need some additional notation. The algebra gl(3) 
is generated by the elements Hi, i = 1, 2, 3 and Eij, = (1, 2), (2, 1), (2, 3), (3, 2), for 
which we use also the notation 



E a — E12, Ep — E23, F a — E 2 i, F/3 — £"32, 

and 

H a — Hi — H 2 , Hp = H 2 — H 3 , H a+f3 = Hi — H 3 . 
They satisfy the relations 

[Hi, Hj] = 0, [Hi, E k i] = (Sik — Su)E k i, [E ai ,F aj ] = 8 ai)0l . H ai , (2.32) 
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E 2 E —OR E E 4- E E 2 — D 

ai^aj. (2.33) 

F 2 F —IE E F A- F F 2 — 

1 Cti 1 Ctj * 1 Cti 1 Ctj 1 Cti T -t ctj 1 oti ~ U ) 

where in the last three relations greek indices and aj take two values a or f3. 
Introduce two-dimensional vectors 

ei = e 2 = (-1, 0), e 3 = ( \, f ), 



(2.34) 



tOl = (1, ^), «% = ( 0,-^), W 3 = (-1,^), 

hat 

e 1 + e 2 + e 3 = 0, + u> 2 + w 3 = 0, (2.35) 



and 

1 2 
e^-^-Wfc); w i = ~^( e j - e k), (2.36) 

where (i,j,k) is a cyclic permutation of (1,2,3). Moreover, 

1 2 
e^j = -a^, WiWj = -a,ij, e i w j = e ij , i,j = 1,2,3, (2.37) 



where a& denotes the (Euclidean) scalar product, a^, i,j = 1,2,3, is the Cartan matrix 
an = 2, aij = -1, i^j, i,j = 1,2,3, (2.38) 



of the algebra 



and 



e ll — ^22 — £33 — 0, 612 — C 2 3 — C 31 — 1, E 2 \ — €32 ~ e 13 — — 1- (2.39) 

It is also useful to define 



£31 = [Fp, F a ], H = w 1 H 1 + w 2 H 2 + w 3 H 3 . (2.40) 

For any ji = (fj,i, // 2 , ^3), such that /i a — n\ — 112 and /ip — 112 — A*3 are non-negative 
integers, denote by 7r M the finite dimensional irreducible representation of gl{3) with the 
highest weight /j, and highest weight vector | 0) defined as 

E a \0)=Ef,\0) = 0, H i \0)= f ji i \0), i = l,2,3. (2.41) 

The dimension of such a representation is dim[7r M ] = \{ji a + l)(/ia + l)(/i a + \ip + 2). 
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Consider the following differential operator 

C = d u - <t>'{u) H - X(E 23 + E 31 + E 12 ) 



(2.42) 



where the components of the vector (f>(u) = [<f>\{u), <p 2 {u)) are the canonical variables with 
Poisson bracket 

1 ,- , x / x u 1 sin(rt-u) 

{M»)Aj(v)} = -S ij e(u-v), e (u) = — + -Y t —^-, 0<u,v<2n ] (2.43) 



n=l 



which are related to U(u) and V(u) by the generalised Miura transform [35 

U(u) = <P' l {uf + <P' 2 {uf + 2<P'[{u), 
y/W(u) = ^0 / 2 ( M ) 3 -20;( M )V 2 ( M )-30 , 1 ( M )0 , 2 / ( M )-0 , 1 / ( M )0 / 2 ( M )-0 / 2 // ( M ). (2.44) 



For the periodic potentials U(u) and V(u), ( |2.18| ) the field <fi(u) is, in general, quasiperi- 
odic 

(f>(u + 2tt) = <t>(u) + 2ixip. (2.45) 

The monodromy matrix associated with the equation 7r M [£]\l/(tt) = (where tt^ denotes 
a representation of gl(3)) is given by 

r-2w 



MJX) = 7T p 



„2nir. 



ipH V exp (A J {V 1 E W + V 2 (u)E 31 + V 3 (u)E 12 )du 

where the symbol V denotes the "path ordered" exponential and 

Vi{u) = e~ 2e ^ {u) , 1 = 1,2,3. 

The associated L-operator 



Lf°(A) 



-wipH' 



M M (A), 



(2.46) 



(2.47) 



(2.48) 



satisfies the "r-matrix" Poisson bracket algebra [40.41 



{LjW (A)®L^(A')} = [v(A/A'),L(f) (A) ® L^(A')] 



(2.49) 



The quantity r w /(A) = 7r M <g> 7iy [ r ] is the "classical r-matrix" given by 

r(A) = -1 4 J2Hi®H t + — -3- V (x-^Eij <g> ^ + A^ <g> ) , (2.50) 

Aa — A 2 — v / 



AS — A" 
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where the second sum is taken over the values (ij) = (12), (23), (31). 
It follows from fl2.49|) that the quantities 



Tf° (A) = Tr WM M / ,(A), (2.51) 
are in involution with respect to the Poisson brackets ( p.l9p ,( p.43|) 



W(A) ) TW(A')} = 0. (2.52) 



Note, in particular, that T^ QQ JX) coincides with Q2.28| ). 



3 Quantum Yang-Baxter equation 
3.1 Universal R- matrix 

The quantum Kac-Moody algebra A = U q (sl(3)) is generated by the elements Xj, jji and 
hi, i = 1,2, 3, and the grading element d, satisfying the commutation relations 

[hi, hj] = , [hi, Xj] = -OijXj , [hi, y 3 \ = a^yj , (3.1) 

q 'hi 

[d,Vi] = Si,3Vi , [d,Xi] = -S it3 Xi , [yi,Xj] = Sij - - , (3.2) 
and the cubic Serre relations 

% J I iQ % J "-'i I J % — 5 

M 3- (3-3) 

yf Vj - [% Vi Vj Vi + Vj vl = °> 

Here the indices i,j take three values i,j = 1, 2, 3; the Cartan matrix a^- is the same as 
(|2.38|) and the "g-numbers" [n] q are defined as 

71 g = — . 3.4 

g-g 1 

The element = h% + /i2 + ^3 is the central element of the algebra A. In the following we 
will always work with the representations of the subalgebra of A with suppressed grading 
element d on which the central element k acts by zero: 

hx + h 2 + h 3 = 0. (3.5) 
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(3.6) 



The algebra A = U q (sl(3)) is a Hopf algebra with the co-multiplication 

8 : A — >A®A, 

defined as 

5(xi) = Xi <g> 1 + q~ hl ® Xi , 

KVi) = Vi®q ht + l®Vi , 
5 (hi) = hi <g> 1 + 1 <g> hi , 

5(d) = d®l + l®c?, 
where i = 1, 2, 3. As usual, we introduce the "twisted" co-multiplication 

8' = ao 8, oo(a®b)=b®a (Vo,6 G A) . (3.7) 

Define also two Borel subalgebras £>_ C A and £>+ C ^4 generated by /i i5 2j, i = 1, 2, 3, and 
/ij, i = 1,2,3, respectively. Let B± be their extensions by d, that is, the subalgebras 
of A, generated by B± and d . 

There exists a unique element [f42| , |4"3| 



ft e B + ® £_ , (3.8) 

satisfying the following relations 

<?'(a) 11 = 11 8(a) (V aei), (3.9) 

(<5®l)ft = ft 13 ft 23 , (3.10) 

(1<8>5) ft = ft 13 ft 12 , (3.11) 

where ft 12 , ft 13 , ft 23 G .4 ® .4 <g> .4. and 1Z 12 = 1Z ® 1, ft 23 = 1 ® ft, ft 13 = (tr ® 1) ft 23 . 
The element ft is called the universal i?-matrix. It satisfies the Yang-Baxter equation 

ft 12 ft 13 ft 23 = ft 23 ft 13 ft 12 , (3.12) 



which is a simple corollary of the definitions ( |3.9| )-( |3TTl| ). The universal i?-matrix is 



understood as a formal series in generators in B + ® £L. Its dependence on the Cartan 
elements can be isolated simple factor 

lZ = q t % (3.13) 
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following where t is an invariant tensor in the tensor square of the extended Cartan 



subalgebra. It can be chosen to be |H, 45] t = hi®h l + h 2 ®h 2 + k®d + d®k 



h 1 = -hi H — ho, h 2 = -ho H — hi, 
3 3' 33' 

and q l satisfies the following defining properties: 

q\yi ® 1) = (jfi <g> q hl )q\ q l (xi <g> 1) = (x< <g> q~ hi )q\ 

(3.14) 

^(l&ij/i) = {q hi ®Ui)q x , q t {l®x i ) = (q~ ht <g> Xi)q l . 

The factor 7?. is the "reduced" universal i?-matrix KcB + ®6_ which is a formal power 
series in yi <S> 1 and 1 <g> x« only 

(with no dependence on the Cartan elements) 

U=\ + {q- q^iyi ®x 1 + y 2 (g)x 2 + y 3 ®x 3 ) + ... (3.15) 

Below we will use the following important property of the reduced universal i?-matrix. 

Proposition 1 Suppose elements I , A{ and Bi, % — 1,2, 3, satisfy the relations 

M = (3.16) 
q-q 1 

Xj = q aij xj A h Bi Xj = q" 0,13 Xj B h (3.17) 
where i, j = 1,2,3. Then 

3 3 

[1®I,TI] = (J2vi®Ai) K-K^yiQBi). (3.18) 

i=i i=i 

The proof of this statement is given in Appendix A. 

3.2 Free field representation and the vertex operators 



The quantum version of the Miura transformation ( 2.44 ) is the free field representation [46 
of the algebra W3, 

-gT(u) = :0 / 1 ( M ) 2 : + :0 / 2 ( M ) 2 :+2(l-^)0 / 1 '( M ) + ^, 
{:j<l)n W(u) = : 0' 2 ( M ) 3 : -3 : #(u) 2 #(ti) : +~G? - 1) : : (3.19) 



2 ' ' ^ ' iV / -rzv / 2 

+ ^O/-l):0i(«)^(«):-|O/-l) a <(«), 
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in terms of a two-component free chiral Bose field <~p{u) = (<pi(u), ^ 2 (w)), 

ip(u) = iX + iPu + \"^ e inu . (3.20) 
t—' n 

The mode operators X = (Xi, X 2 ), P = (Pi, P2) and b n = (&i jn , 62,™) (n = ±1, ±2, . . .) 
satisfy the commutation relations of the Heisenberg algebra 

[Xi,Pj] = -^Sij, [bi,n, Km] = -Y$n+m,0 Sij, i,j = 1,2 (3.21) 

where the parameter g is related to the central charge c as 

50 - c - J (2 - c)(98 - c) , . 

9 = 4g "> c = 50-2A(g + g- 1 ). (3.22) 



Let J-p be the highest weight module over the algebra (|3.21|) , with highest weight vector 
I p)i P — (P11P2) defined as 

Pi\p)=Pi\p), a i<n \p) = 0, n>0, i = l,2. (3.23) 

For generic c and p, the Fock space T v is isomorphic to the space Va 2 ,a 3 with the highest 
weights 

a P1+P2 1 c-2 2p 2 (g| - 3gf) 



The space 



•^"p © p+g(niei+n2e2)) (3.25) 

ni ,ri2=— oo 



supports the action of the vertex operators 

= g5 : e -2e^H . ? « = 1,2,3, (3.26) 

defined as 

oo , OO , 

. e -2 e M u) .= exp ( _ 2 y exp / _ 2le ^ x + p M)) exp ( + 2 y ^ e -i»«) 

n=l n=l 

(3.27) 



where the root vectors are given in ( 2.34j) and 



q = e iwg . (3.28) 

Using ( |3.21| ) it is easy to check that 

V l {u 1 )V j {u 2 ) = q a ^Vj(u 2 )Vi(ui), Mi > u 2 , 
[P,V»] = -^-V-H, 
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(3.29) 



where the ay are given by (|2.38 ) and 

Vj{u + 2vr) = q- 2 e-^ p ^V 3 [u). 
It will sometimes be convenient to use the operators 



Zj = e 



ziz 2 z 3 = l, j = 1,2,3, 



(3.30) 



(3.31) 



where the vectors Wj are defined in ( 2.34| ). For example, with these operators the relations 
(|3T29|),(|3T30|) take the form 

ZiVj{u) = q 2e ^V j (u)z i , (3.32) 



where is defined in ( j2.39| ) and 



Vi(u + 2ir) = q 2 zj 1 z k Vi{u), 



where k) is a cyclic permutation of (1, 2, 3). 



(3.33) 



3.3 The universal L-operator 

Consider the following operator f47|,f48 



2- 



where 



£ = e l7T ™Vexp ( j K(u)du 



h = -(e 1 h 1 + e 2 h 2 + e 3 h 3 ) 



with the vectors e\ y e 2 and e 3 defined in ( p.34| ) and 

K(u) = V 1 {u)y 1 + V 2 (u)y 2 + V 3 (u)y 3 . 



(3.34) 



(3.35) 



(3.36) 



Here hi, y iy i = 1,2, 3, hi + h 2 + h 3 = 0, are the generators of the Borel subalgebra 
£>+ C U q (sl(3)) defined above. The ordered exponential in ( |3.34| ) (the symbol V denotes 
the path ordering) is defined as a series 



C = e inPh Y^ / K(u l )K(u 2 )---K(u n )du l du 2 ...du n . 

n = 27r>Ul>M 2 >-">«n>0 



(3.37) 



The integrals in ( |3.37|) converge when 



oo < c < -2, 



(3.38) 
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(which corresponds to the range < g < 2/3). For the values of c outside this quasiclas- 
sical domain, ( 3.38 ), the integrals in ( |3.37|) can be defined via analytic continuation in 
c 0. The operator ( 3.34)) is an element of the algebra B + whose coefficients are operators 
acting in the quantum space (|3.25| ). Following the arguments of j| it is not difficult to 
show that this L-operator satisfies the Yang-Baxter equation 



C)(C®1)K. 



(3.39) 



where 1Z is the universal .R-matrix for the algebra U q (sl(3)) defined above. The only 
property of the universal R- matrix required for ( |3.39| ) to hold is relation ( |3.9| ). 

An important remark is in order. As the reader might expect, the L-operator ( |3.34|) 
can be obtained as a specialization of the universal .R-matrix and the equation ( |3.39| ) 
is, in fact, a special case of the Yang-Baxter equation (|3.12| ). The exact correspondence 
(first established in |I| for the case of U q (sl(2))) is formulated as follows. Each side of the 
equation ( |3.12| ) is an element of the direct product B + <8> A <8> £>-• Therefore to reduce 
(|3.12|) to (|3.39| ) it is enough to find an appropriate realization of the Borel subalgebra 
B- (appearing in the third factor of the product) in the quantum space (|3.25| ). Let us 
identify the generators Xi G B-, i = 1, 2, 3, of this Borel subalgebra with the integrals of 
the vertex operators 



q-q 



-i 



2- 



Vi(u)du, 



1,2,3. 



(3.40) 



One can check that these generators satisfy J49[ the Serre relations (|3.3j ). Substituting 
(|3.40|) into the reduced universal .R-matrix, ( |3.15|) , one gets an element of B + whose 
co-ordinates are operators acting in the quantum space ( |3.25| ). As expected, it exactly 
coincides with the "P-ordered exponent from ( |3.34|) . 



Proposition 2 The specialization of the "reduced" universal R-matrix 1Z G B + ® i3_ for 
the case when X{ G B-, i — 1,2, 3, are realized as in ( \3.4C\ ) coincide with the V -exponent 



n 

Xi=(3.4Q) 

where lC{u) is defined in ( \3. 3(Jj ). 



V exp 



2t 



fC{u) du 



(3.41) 



This is precisely the conjecture stated in ref. 0]. It was formulated there for U g (sl(2)), 
but the statement is rather general and (with obvious modifications) holds for any quan- 
tized Kac-Moody algebra. A proof of the conjecture of @] based on combinatorial prop- 
erties of the representations of quantum groups has been obtained by E.Frenkel and 
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F.A.Smirnov |50|. Here we present an alternative (and, in fact, rather elementary) proof 
of ( 3.41 ) which is based on the defining properties of the universal i?-matrix. 

The reduced universal i?-matrix ( 3.15|) is a formal power series in the generators 
Vi — Vj ® 1 an d Xi = 1 <S) Xi, i = 1, 2, 3, such that xjy~j = TjjXi, therefore we can write it as 
1Z = IZijj^Xi). Then, substituting (|3.13|) into ( |3.11|) one obtains 

TZiyj, x\ + x'l) = IZiy^x't) Tl^^x'l) 



where y i — yi ® 1 Cg) 1, x[ — 1 ® q h% ® x^ x" = 
Serre relations ( |3.3|) of £L , commute with y • 



(3.42) 

1 <g) Xi ® 1. Both sets x\ and x'[ satisfy the 



x iUj UjXii x^yj yjXi 



(3.43) 



and also obey the relations 

x\x] = q a ^x]x\ (3.44) 

Conversely, we can regard ( |3.42| ) as an identity for IZiij^ Xi) which should be valid once all 
the arguments involved there are subjected to the above relations (namely, (|3~3|), Q3.43D 
and (pi). 

Consider integrals more general than ( p.40| ) which also satisfy the Serre relations ( |3.3|) 

1 f Ul 

Xi = -r / Vi(u)du, i = 1, 2, 3, (3.45) 

9 - 9 Ju 2 

where 2tt > u\ > u 2 > 0. Denote by C(u\, u 2 ) the result of substituting these generators 
Xi G £L into the reduced universal i?-matrix 



£(wi,w 2 ) = K{y h Xi 

Xi=(3A5) 

Note that from ( |3.15| ) it follows that for (u\ — u 2 ) — > 

/-Ml 

£(u u u 2 ) = l + K(u)du + 0((u 1 -u 2 ) 2 ) 



(3.46) 



(3.47) 



«2 



where K(u) is defined in 

Now, let 27r > u\ > u 2 > M3 > 0, then the identity ( [3.42| ) implies that 

C{u\, u 3 ) = C(u 1} u 2 )C(u 2 , u 3 ), 

since two sets of quantities 

1 f Ul 
Xi = -r I Vi(u)du, 



(3.48) 



9-9" 



"2 
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x" = — - — - / Vi(u)du, 
q — q~ l 



satisfy all the required relations ( |3.3| ), ( |3.43| ) and ( |3.44| ) . In particular, the relation ( |3.44| ) 



readily follows from ( p.29|) . With an account of (|3.47|) the functional equation ( |3.48|) has 
a unique solution 

C( Ul , u 2 ) = Vexp(^J K{u) du^j . (3.49) 
Finally, setting u 2 = and u\ = 2ir one arrives at the statement of Proposition 0. 



4 Quantum T- and Q-operators 
4.1 Symmetry relations 

Consider the symmetry properties of the L-operator ( |3.34| ). Any linear orthogonal op- 
erator U in two-dimensional Euclidean space induces an automorphism U of the two- 
component Heisenberg algebra ( [3.21|) by the rule U[ip(u)] = XJip(u). In particular, the 
Weyl reflections associated with the vectors Wi, i = 1, 2, 3, generate the symmetric group 
S3 permuting the vertex operators 

VijWiiu)] = Vj(u), aiAVj{u)\ = Vi(u), <JiAVk{u)] = V k (u), (4.1) 

<Tij[zi] = zj 1 , aijizj] = z~\ Oij\zk\ = z k \ (4.2) 

where (i,j,k) is any permutation of (1,2,3). Define also the cyclic automorphism r = 
0"i20"23 corresponding to the rotation of the Dynkin diagram of 

r[V t (u)} = V l+1 {mod3) {u). (4.3) 

The weights Q3.24 ) can be written in a symmetric form 



. r x r 2 + r x r z + r 2 r 3 c-2 r x r 2 r 3 
A2 " g + ^4~' As " ' (4 ' 4) 

where 

r { = WiP, i = l,2,3, (4.5) 
with Wi given in ( 2.34j) . From the symmetry properties 



T[ri] = r i+ i ( m od 3), o\?N = -rj, °iA r A = a iA r k] = ~ r k, (4.6) 
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it is obvious that the cyclic automorphism r transforms the Fock module T v to an iso- 
morphic module 

r\T v \ ~ T v , (4.7) 
while the Weyl reflections negate p and the weight A 3 

criA^Fp] ~ T- v ~ Va 2 ,-a 3 - (4.8) 

Further, the algebra of outer automorphisms of B + C U q (sl(3)) also contains the 
symmetric group S3: <j[yi] = Va(i), &[hi\ = h a (i)- Again by r we denote the automorphism 
corresponding to rotation of the Dynkin diagram of 

T [yi] = Vi+l (mod 3), r[hi] = h i+1 ( mo d 3) , (4.9) 
and by the transpositions 



VijiVi] = Vj, VijlVj] = Vu Vijlyk] = Vk, 

aij[hi] = hj, CiAVh] = ^, Oij[h k ] = h k . 



(4.10) 



The L-operator (|3.34|) is invariant with respect to the diagonal action of the automor- 
phisms 

(a®a)[L\ =£, a G S 3 . (4.11) 

4.2 The T-operators 

We will need several specializations of the universal L-operator (|3.34|) corresponding to 
various matrix representations of the Borel subalgebra B + C U q (sl(3)). First consider the 
so-called evaluation maps from £>+ to the finite dimensional algebra U q (gl(3)). This latter 
algebra is generated by the elements Hi, i = 1, 2, 3 and Eij, = (1,2), (2,1), (2,3), (3,2), 
for which we use also the notations 

E a = E12, Ep = E 23 , F a = E 2 i, Fp = E 32 , 

and 

H a = Hi — H 2 , Hp = H 2 — H 3 , H a+ p — H\ — H 3 . 
They satisfy the relations 



[Hi, Hj] = 0, [Hi, E kl ] = (5 ik — 5u)E, 



kl ■ 
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^ - [2]g -Ea, + E aj E^. = 0, 

where = (q n — q~ n ) / {q — q^ 1 ) and the greek indices and in the last three relations 
take two values a or /?. The algebra U q (gl(3)) has an outer automorphism CX13, induced 
by the automorphism of its Dynkin diagram: 



0"i3[-Ea] = Ep, ^n[Ep] = E a , cr 13 [F a ] = Fp, cmlFp] = F a , 

(4.13) 

<7 13 [#i] = ~H 3 , a 13 [H 2 ] = -H 2 , a 13 [H 3 ] = -H x . 
Let us introduce the following root vectors: 

F a+ p = q^FpF a - q~^F a Fp, F a+ p = q'^FpF a - q^F a Fp, 

and 

£ 3 i = F a+ pq H ' +H \ E 31 = F a+ pq~ H ^ H \ 
There are two non-equivalent evaluation maps B + — > U q (gl(3)). The first one is 

£v t (y 1 ) = E 23 , £v t (y 2 )=tE 31 , £v t {y 3 ) = E 12 , (4.14) 

£v t (h 1 ) = H 2 -H 3 , £v t (h 2 ) = H 3 -H x , £v t (h 3 ) = H x - H 2 . (4.15) 

where t is a spectral parameter. The second evaluation map £v t : B + — > U q (gl(3))[t] has 
the form 

£v t = 0- 13 • £v- t ■ 0"13, (4.16) 

where the automorphisms a\ 3 entering this formula are defined in ( |4 . 1 0[ ) , (|4.13|) . Explicitly 

one has 

£v t (yi) = E 23 , £v t {y2) = tE 31 , £v t (y 3 ) = E 12 , ( 4 -17) 

£^ t (h 1 ) = H 2 -H 3 , £^ t (h 2 ) = H 3 - H u £^ t {h 3 )=H x -H 2 . (4.18) 

One can easily check that both these maps are algebra homomorphisms as all the defining 
relations (|3.3[), ( |3.1| ) of the algebra £>+ become simple corollaries of ( |4.12| ). 

For any fi = (fii, fi 2 , fi 3 ), such that n a — ^\ — \i 2 and \ip = \i 2 — \i 3 are nonnegative 
integers denote by 7r M the finite dimensional representation of U q (gl(3)) with the highest 
weight /i and highest weight vector | 0) defined as 

E 12 \0) =E 23 \0) = 0, H i \0) = f JL i \Q), 1=1,2,3. (4.19) 
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Further let ir^t) be the representation of B + given by the composition of 7r M with the 
evaluation homomorphism £v t . Evidently, the operator- valued matrices 



L fl {t)=n ll (t)[£], 



(4.20) 



(whose matrix elements are operators acting in (|3.25|) ) satisfy suitable specializations of 
the quantum Yang-Baxter equation (|3.39| ). These operators (|4.20|) are quantum analogues 



of the classical L-operators (|2.48|) whereas the operators 



T^t) = Tr 



4 



(4.21) 



are the quantum analogues of T^(A) in (|2.51| ) (the spectral parameter t is related to A 
in ( |2.48 ) and (2.51) as t — A 3 ). From the standard arguments based on the Yang-Baxter 
equation it follows that 

[T,(t),TV(t')] = 0. (4.22) 
Let us show that they also commute with all local IM (|2.8|) , 



[TV(t),I fc ] = 0. 



(4.23) 



The characteristic property of the local IM is that their commutators with vertex operators 
(|3.26|) reduce to total derivatives |54 



[lk,Vi(u)] = d u{ ■ 0®(u)Vi(u) : } = duQfiu), 



(4.24) 



where 0^\u) are some polynomials with respect to the field d u <f>(u) and its derivatives. 
Obviously, the operators Q^\u) in Q4.24Q obey the same periodicity properties as the 
vertex operators Vi(u), ( |3.30| ), 



G%\u + 2vr) = q- 2 e- 4inP ^Q^\u) 
and satisfy the commutation relations 



(4.25) 



U\ > u 2 , 



-g e,- 6 



(i), 



(4.26) 



where ay is defined in ( |2.38|) . It is easy to check that the triple 



(0/ 



(4.27) 
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satisfies the conditions (|3.16|) , (|3TT7|) of Proposition [l| with the generators Xi G £>+ realized 
as in (|3.40 ) Therefore taking into account ( 3.41 ) and applying ( 3.1 8|) one obtains 



[I k ,C} = Q k (2n)C-Ce k (0), (4.28) 

where 

3 

e fc (tO = X>e<°(«)- (4.29) 



i=l 



Now substitute ( |4.21| ) into the RHS of ( ^4.23| ). Using the relation ( |4.28| ), the commutation 



relations fl3.1|) , ( |4.26| ) and the cyclic property of the trace it is not difficult to show that 
the commutator (|4.23| ) actually vanishes. 



The most important T-operators correspond to the three-dimensional (fundamen- 
tal) representations of 7t Wq , and 7r W(3 of U q (gl(3)) with the highest weights oj a = (1,0,0) 
and up = (1,1,0) (The corresponding representation matrices are given explicitly in 
Appendix |B|). We will introduce special notation for these T-operators (adjusting the 
normalisation of the spectral parameter for later convenience) 

T(t) = T ( i,o,o)(*g-5), T(t) ee T (1 , 1>0 )(*g-3). (4.30) 

Computing the trace in ( |4.21| ) in this case one obtains 

oo oo 

T(t) = £fG ni T(t) = tn G n - (4.31) 

n=0 n=0 

Here 

G = z 1 + z 2 + z 3 , G = z{ x + z 2 l + ^ , (4.32) 
with zi, z 2 and z 3 defined in ( |3.31|) and for n > 1 



r — y 4- r 7"^ 4- r ^ 

— ^1 ^312 ~r ^2 ^123 ~r ^3 J 231) 

1 NnfT -1 r(n) -1 r(n) -1 x(n) 

"J-; W — ^213 ^2 ^321 "r ^3 J 132 



(4.33) 



4* = / ^(^il^^lVihl^W^W^W- (4.34) 

■ ■ ■ Vi(u 3n ^2)Vj(u 3n ^i)V k (u 3n ) J dlXi d« 2 • • ■ ^3n, 

where (i, j, k) is a permutation of (1, 2, 3). 
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It follows from ( |4.22j ) and (|4.23|) that the operators G n and G n can be regarded as 
the "nonlocal" IM which commute among themselves and with all the local IM ( |2.8|) 



[G fc , Gj] = [G k , Gi] = [G k , G{\ = [G k , I,] = [G k , h] = 0. 



(4.35) 



Definitions similar to those in ( 4.21 ) but with the map ( |4.18| ) instead of ( |4.15| ) lead 
to a new set of T-operators 



T„(f) = Tr= it) [e"* h JC\, 



(4.36) 



where 7F^(t) denotes the representation of B + obtained by the composition of the represen- 
tation 7r M of U q (gl(3)) with the evaluation map fl4.|18| ). The operators T /t (t) with different 
values of t commute among themselves and with all T M (£), 



[T,(t),T, l (f)] = [T,(t),T, / (t')]=0. 



(4.37) 



In general for the same 7r M the formulae ( |4.21| ) and ( [4.36|) define different operators. 
However, when 7r M is restricted to symmetric powers of a fundamental representation -K ula 
and 7r^„ then the two sets of T-operators coincide (up to a rescaling of spectral parameter). 
For future references it will be convenient to define 



■ (m,0,0) 



(tq- 



(4.38) 



and 



TT m (t) = T( m>mi0 ) (tc? m+2 , 



T( m , m ,o)(^ m " 5 ), (4.39) 

where m = 0, 1, 2, . . . , oo. In particular, T (t) = T (t) = 1 while Ti(t) = T(t), Ti(t) = 
T(t) are exactly the same as in ( [4.31]) . 



4.3 The Q-operators 

The operators (|4.21|) and (|4.36|) are obviously the CFT versions of Baxter's transfer- 
matrices of the lattice models related to the U q (sl(3)) algebra. It is well known that 
the theory of solvable lattice models involves other important objects, namely Baxter's 
Q-matrix. Similarly to the standard CFT case appropriate W3 CFT versions of Q- 
operators can be obtained as certain specializations of the universal L-operator ( |3.34|) . 

The generators jji, hi, i = 1,2,3, hi + h 2 + h 3 = 0, of the Borel subalgebra B + C 
U q (sl(3)) satisfy the commutation relations (|3.1|) 

[hi, 11,1 = 0, [hi,yi\ = a i5 , (4.40) 
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and the Serre relations ( |3.3| ). The later can be written as 

ViVij = Q^VijVi, h3 = 1, 2, 3, % ^ j, (4.41) 
where j/jj stands for the g-commutator of yi and yf 

Vij = ViVj - IViVi- ( 4 - 42 ) 
Let &f ,sl~ ,Hi, i — 1,2 be the generators of two independent q-oscillator algebras H g , 

[Hi, at] = ±af, qa+ar - q^ar&f = . (4.43) 

Define the following homomorphism p t : B + — > Y\ q ® Y\ q 

Pt{Vx) = Q*t(q- g -1 ) g~ Wa a7/aJ, p t (y 2 ) = a 2 , p t (y 3 ) = afq n2 , 

(4.44) 

Pt(h x ) = -Hi + H 2 , Pt (h 2 ) = -H 1 -2H 2 , Pt(h) = 2H 1 + H 2 , 

where t is a spectral parameter. Note that this map subjects the generators ?/j, i — 1, 2, 3 
to some additional relations which are not implied by (|4.40| ) and the Serre relations Q4.41]) . 

i 

q 2 t 

Pt (qymi - g~ 1/212/3) = pt (W32/2 - q~ 2/22/13) = 7, pt (2/23) = 0, (4.45) 



where y^ is defined in ( |4.42|) . Define also the maps pi(t) and Pi(t) : B + — > \-\ q <g> H 
i — 1,2,3 obtained as compositions of ( |4.44j ) with the automorphisms ( (4.9| ),( 4~TU ) of the 
algebra B + 



Pi(t)= Pt-r~ l+ \ p^t) = p- t ■ r- i+1 ■ a jk 2 = 1,2,3, (4.46) 



where pt is given by (|4.44 ) and p^t is obtained from ( |4.44| ) by the replacement t — > —t. 
Let 7Tj and 7fj be any representations of H g ® H g such that the traces 

Z t = Tr 7TiPi{t) [e 2 ^ h ] and Z t = Tr— pM [e 2 ^% (4.47) 

converge and do not vanish for 

3m e x p < 0, 3m e 3 p < 0, (4.48) 

where p = (pi,p 2 ) G C 2 . 
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Define the operators 



r-l 



and Ai(t) = Z i 1 Tr__+ (t) [e* 



[e lnPh C], (4.49) 



where C is given in (|3.34|) . Since we are interested in the action of these operators in 
T v the operator P in (|4.49|) can be substituted by its eigenvalue p = (pi,P2)- The 
definitions ( |4.49| ) obviously apply to the domain ( |4.48| ). However they can be extended 
for all complex p e C 2 (except for some set of singular points for real values p\ and P2) 
by an analytic continuation in p. 

Obviously, the operators ( |4.49| ) can be written as a power series. For example, 



A,(«) 



E E • 

n=l {cr,=l,2,3} 



42< 



(Ti, a 2 , ■ ■ ■ , cr 3n )J 3n ((Ti, cr 2 , . . . , cr 3n ), 



where 



and 



a 3n (0-1,0-2, 



0- 3r 



Kx («l) (w 2 ) • • • V a3n (u 3r , 



o 3 „. 



Zr'Tr, 



"I/a 



(4.50) 



(4.51) 



(4.52) 



Note that these coefficients vanish unless the "total charge" of the product of vertex 
operators in ( [4.51| ) 



,(0 



(4.53) 

0(t n ) so that (POP 



is zero. The spectral parameter t is absorbed in the coefficients a\ n 
is a power series in t. 

A remarkable feature of the definitions ( [4.49 ) is that the coefficients ([4.52 ) are com- 
pletely determined by the commutation relations ( [4.40|) , ft4.41|) , (|4.45| ) and the cyclic 
property of the trace, so the specific choice of the representations 7Tj and Tfj is not signifi- 
cant as long as the above convergence property is maintained. The operator coefficients 
in the power series (|4.50| ) can be expressed in terms of basic nonlocal IM (|4.33| ). For 
example, calculating the first nontrivial coefficient in ( |4.50| ) one gets 



Mt) 



Mt) 



q2(qz i G 1 - Gi 



(q 2 - l)(zj - q 2 Zi)(z k - q 2 Zi) 



q^Zid - gGi 



(q 2 - l)(q 2 Zj - Zi)(q 2 z k - z { ) 



t + 0{t 2 ), 



t + 0(t 2 



(4.54) 



27 



where Zi, i = 1,2,3, defined in (|3.31|) . For further references introduce an alternative 
set of nonlocal IM, H t (n) , H; , i = 1,2,3, n = 1, 2 . . . oo defined as coefficients in the 



expansions 



where 



log Aj(t) = -J2 s n , A, (t) = - s n 



(4.55) 



n=l 



n=l 



s = t(g~ 1 T(l-g)ft. (4.56) 

The operators Aj(t) and Aj(£) simplify considerably when acting in the space T v with 
values of 2p belonging to the weight lattice of the algebra sl 3 



2p = n\Wi + n 2 w 2 , ni,ri2 G Z. 
In this case the coefficients ( |4.52 ) can be written as 

a3n(cri, 02, ■ ■ ■ , 0"3n) = Tl u (u ai U c 

where u\, u 2 are generators of the Weyl algebra 

uiu 2 = qu 2 u 1 , 



with the trace defined as 51 



and 



Tr u «u™) = 5 niQ 5 m ,o, n, m e Z 
i 

qz t 



u 3 



(q-q x ) 



1X3 U 2 V 



(4.57) 



(4.58) 



(4.59) 



(4.60) 



(4.61) 



Obviously, the expression (4.58) is invariant with respect to the cyclic automorphism r 
in ( fO| ) therefore all three operators Aj(t), % = 1,2,3 coincide. With the help of ( |4.5(J| ), 
(|4.58|) they can be written as 

r-2TT 



2p=(4.57) 



Tr 



Pexp 



u 1 V 1 {v)+u 2 V 2 {v)+u 3 V 3 (v)yv\, t = l,2,3. (4.62) 



Exactly the same expression but with u\ and u 2 interchanged (which is equivalent to the 
replacement q — > q' 1 in fl4.59|) ) holds for the operators Aj(t). 

The Q-operators (the CFT analogues of Baxter's Q- matrix for U q (sl(3) related lattice 
models) are defined as 



Qi(t) = t^r Ai(t), Qi(t) = t ir Ai(t), 



(4.63) 
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where Wi, i = 1, 2, 3, are defined in (|2.34|) . 

Evidently, the operators Qi(t) and Q, ; (t) with different spectral parameters commute 
among themselves and with all the operators T M (t) and T^(t) constructed above. They 
also commute with all the local IM defined in (|2l 



[Qi(*),I fc ] = 0, [Q,(t),I fc ] =0. (4.64) 



Indeed, the proof of the similar statement ( |4.23| ) for the T-operators given above does 
not depend on a particular choice of the representation of B + and relies only on the 
commutation relations (|4.40|) , (^4.41|) and the cyclic property of the trace. Therefore this 
proof applies to ( |4.64| ) as well. 

As is known |52]-|54| the local IM do not change under the transformation g — > g^ 1 , 



provided one makes simultaneous substitution 4>{u) — > g 1 4>{u) 

l k {g,4>{u)} = l k {g-\g- l cj>{u)}. (4.65) 

Obviously, the nonlocal IM G n and G„, (as well as H- n ^ and H^) do change and there 
exists an infinite set of dual nonlocal IM Q obtained from (|2.8| ) by the above substitution. 
For instance, 

F< — ~3 7(") I ~3 J( H ) I r 9 7( n ) (A ««^| 

u « — z l ^312 ' z 2 J 123 "r z 3 J 231' I^.UUj 



J(n) 
J ijk 



(ViMVjMVkfo) Vtiu^VjMVkiue) ■ ■ ■ (4.67) 



^(M3n-2)^(M3n-l)Vfc(M 3n ) ) du X du 2 . . . du 3n , 



where k) is a permutation of (1, 2, 3) and 

Viiu) =qh- 2e ^ {u)/9 , q = e in/g . (4.68) 

Similarly, define 

%{g, 4>{u)} = G n {g-\ g-^iu)}, (4.69) 

and 

^{g,^)} = Hf^- 1 ,^ 1 ^)}, H^lg, 4>(u)} = H- n) {g _1 , g~ l <f>(u)}. (4.70) 
For g < 3/2 the above definitions require analytic continuation in g as described in [^J. 
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4.4 Vacuum eigenvalues 

The Fock space ( |3.25|) naturally splits into level subspaces IFp\ I 

oo 



0,1,2, 



defined as 
(4.71) 



All the IM (both local and nonlocal) act invariantly in each level subspace Tp\ In 
particular the vacuum state is an eigenstate of all these operators, 



Ife \p) 



[(vac) 



P), 



G n \p) = G^\ P ), 



Gn\P)=G n \p). 



Below we list a few vacuum eigenvalues jj^ ac ^ which readily follow from ( 2.16Q 



j(vac) 


= A 2 - 


c 

24 ' 




1 2 


= A 3 , 






j(vac) 


= (A 2 - 


c + 6, 
24 ' 


A 3 , 


j (vac) 
1 5 


= A^ + 


^- 

3 3 


C + 8 Aj 
8 2 



(c + 15) (c + 2) c(c + 23) (7c + 30) 
— A 2 



192 



96768 



(4.72) 



(4.73) 



Using ( |4.33| ) one can express the vacuum eigenvalues of the nonlocal IM as Coulomb- 
type integrals in a circle. In general they cannot be evaluated in a closed form. Even the 
simplest eigenvalues G^ ac ^ and G^ reduce to the generalised hypergeometric integral 



J{x,y) 



1 dtoi f 1 dw 



. x+g ,y+g 



U)\ Jo UJ 2 (1 - ^l) 9 (l - ^2) 3 (1 - UiU 2 ) 9 

T(g + x)T(g + y)T 2 (l~g) f g, g + x, g + y _ 
T(l + x)T(l + y) 3 2 V l + x, l + y ' 



(4.74) 



For example, 



(vac) 



Z2 



2- 



du?, 



«3 + 27T PUl 

du\ / du 2 - 



-2ip(eiu\ + e2U2+esus) 



(2 sin 



Ui - u 2 



)(2sin 



)(2sin 



u 2 -u 3 , 



(4.75) 



can be expressed through (|4.74j) in the following alternative forms 

q(voc) 



3 

= — 2nq 2 


g(g 2 zi 


-Z2, 


J32 


- (Q 2 Z2 


-z 3 ] 


Jl2 


3 

= — 27rg 2 


q(q 2 z2 


~z 3 , 


Jl3 


- (<7% 




J23 


3 

= -2nq 2 




~ zi] 


J21 


- (<l 2 zi 


- Z2, 


J31 



(4.76) 
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where 

J i3 = J(2e iPl -2e jP ), (4.77) 

with J(x,y) defined in (|4.74|) and Z\, z 2 and z 3 denote the vacuum eigenvalues of the 
corresponding operators in ( |3.31 ). The vacuum eigenvalue of Gy ac is obtained from 



q{vo,c) ^ ^ e negation of the overall sign in ( 4.75Q , ( |4.76| ) and the replacement 



z i — > z 3 j z 2 z 2 ) 

induced by the automorphism oi 3 , 



Z 3 -> 



ei e 3 , 



e2 -> e 2 , 



(4.78) 



(vac) 



= 2nq~ 


3 
2 




- -j 1 : 


Jl2 


- (?V 


-^r 1 


)j 32 


= 2nq- 


3 
2 




- -r 1 : 


Jsi 


- (sv 






= 2nq~ 


3 
2 


q(q 2 z^ 


- ^: 


J23 


- (?V 




)Jl3 



(4.79) 
(4.80) 



Using (|4.54j ), ( |4.76|) and ( |4.79|) one can obtain the vacuum eigenvalues of the nonlocal 
IM H 4 (1) and uf } defined in flOg|) , 



(i)(i>ac) 



7T# 3 J, 



T 3 (l - 0) sin(7r#) 



"7rW(«ac) 



7T # 3 Jjfc 



T 3 (l - g) sm(ng) 



(4.81) 



where (z, j, k) is a cyclic permutation of (1, 2, 3). 

It will be useful to consider a special choice of the vacuum parameter p = (p, 0). In 
this case, we have 



V 
2 



e 2 p) = -p, 



e-iP) 



p 
2' 



and 



Zx = z- 1 = e 27rip , 



z 2 = l. 

Both G < f ac ' 1 and H^ <yVac ' 1 simplify considerably in this case, 



G { r c \p, 0) 



47r 3 r 2 (i-^)r(i-f 



T(l - f )r(l -g + p)T(l -g- p)T(l - g/2 + p)T(l - g/2 - p) ' 



(l)(vac) 



(p>o) 



7T 



<? 3 r (1 - f ) 



T(g+p)T(g/2 + p) 



2 sinM r (1 - f ) T (1 - g) r(l - g + p)T(l - g/2 + p) ' 
where we have used eq. (4.4(5)) of to simplify the hypergeometric functions J^. 



(4.82) 
(4.83) 

(4.84) 
(4.85) 
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5 Functional Relations 



The T- and Q-operators introduced in the previous Section enjoy a variety of important 
functional relations (FR). In fact, many such functional relations have been previously 
obtained in the context of the solvable lattice models associated with U q (sl(3)) algebra 
(see, e.g. fll2]-|l4], |16|, [L7| ) . It turns out that the explicit construction of the Q-operators 
given in the previous Section, (which defines them as special quantum transfer matrices 
associated with the g-oscillator algebra), can be used to considerably simplify the deriva- 
tion of these FR. In fact, we found it rather remarkable that all known FR relevant here 
(as well as some new ones) can be obtained as elementary corollaries of just two basic 
relations which express the operators T At (t), ( |4.21| ), and T M (t), ( [4.36| ), corresponding to 
finite dimensional representations tt^ of U q (gl(3)) (i.e. those with integral non-negative 
weights n a = n\ — fi2 > and /ig = \ii — Hs > 0) in terms of the operators Qi(t) and 



z 1 det 



z Q det 



Qi(tq 



where fx\ — \i\ + 1, /i' 2 = /i2, h'z = ^3 — 1, an d 

zq = {zi - z 2 )(z 1 - z 3 )(z 2 - z 3 ), 



i,i=l,2,3 



(5.1) 



(5.2) 



(5.3) 



with zi, Z2 and z 3 defined in ( |3.31|) . The operators Qi(t), Qj(t) are not functionally 
independent. They satisfy the relations 



CiQM = Qj{tq)Qk(tq- L ) - QkWQjitq- 1 ), (5.4) 
-CiQM = QjWQ&q-^-QkWQjitq- 1 ), (5.5) 

where (i,j,k) = cycle(l, 2, 3), and 

i 

Ci = z?(zj - z k ). (5.6) 

The proof of these relations is given in Appendix B (see also the discussion at the end 
of this Section). Below we present some simple corollaries of ( |5.1|) , ( |5.2|) , ( |5.4j) , and (|5.5| ). 



The operators T(t) and T(t) defined in ( 4.31 ) corresponding to the three-dimensional 
fundamental representations of U q (gl(3)) satisfy the following relations with the Q-operators 

Q^-Titq^Q^+Titq-^Q^tq-^-Q^tq- 3 ) = 0, (5.7) 
Q^-T^^Q.^+T^-I)^.^- 1 )-^^- 3 ) = 0, (5.8) 
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where i = 1,2,3, These relations are obviously the ?7 g (s/(3))-analogues |1| of Baxter's 
famous T-Q relation [11]. Given T(t) and T(t) each of the equations Q5.71 ) is a third 



order difference equation for Q's which should have three linearly independent solutions. 
Since T(t) and T(t) are single- valued functions of t, (i.e. periodic functions of logt) the 
operators Qj(t) and Qj(t) can be interpreted as the "Bloch wave" solutions of ( W^t ) which 
satisfy the "quantum Wronskian" condition 



zq = det ||Qi(t<7" 



-2^ 



Uj'=l,2,3 



det ||Q,(tg : 



2.n 



l*,i=l,2,3 



(5.9) 



which is just a particular case of ( |5.1| ), ( |5.2| ) for /Xj = 0. 

The equations (|5.7|) and (|5.8|) imply more familiar expressions 



T(t) 



Q(tq*) , Q(tq->)Q(tq*) , Q(tq 



Q(tq> 



+ 



T(t) = gW 



Q(tg2)Q(t g - 
Q(tg-f)Q(tg 



+ 



Q(*?" 
Q(tg- 



(5.10) 



(5.11) 



Q(tga) Q(tg2)Q(tg- 5 ) Q(*grj) 

where {Q(t), Q(i)} is any of the six pairs {Qj(t), Qj(£)} with z 7^ j. More general expres- 
sions Jl6| of this type for the T-operators with arbitrary weight fi can be obtained, for 
example, from ( |5 . 1 U| ) , ( [5.1 1| ) and the relation ( p.lfiQ given below. 

Further, using ( |5.4| ) and ( |5.5| ) in fl5.1| ) and ( p.2| ) one can easily express the operators 
(1P8D and (I^D as 

3 

T m (t) = z l ^ Q Q l (tg m+ i)Q 4 (tg- m ^), 



T m (t) 



8=1 

3 



(5.12) 
(5.13) 



8=1 



(5.14) 
(5.15) 



Using these two results, the well known fusion relations, \Tf\, can be deduced, 

T m (tq)T m (tq- 1 ) = T m (t) + T m ^(t)T m+1 (t), 

T m (tq)T m (tq- 1 ) = T m (t) + T m -i(t)T m+1 {t), 

Note that the formulae ( |5.1| ) and ( [5.2| ) strongly resemble the first Weyl formula for 
characters of gl(3). Note also, that the corresponding analogue of the second Weyl formula 
for the transfer matrices (sometimes called the quantum Jacobi-Trudi identity ]T7|]) has 
been known for some time [16| 



■ (Ati,M2,0) 



(t) = det 



(5.16) 
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where 

T (°)(t) = T (3)( t ) = i, T (D(t) = T (1A0 )(t), r^t) = T (1 , 1>0) (t), (5.17) 
and T^ k \t) = 0, for A; < and /c > 3. The integers /if, i = 1, . . . , /ii, read 
/if = 2, 1 < z < [M 2 ; /if = 1, /i 2 < « < /ii- 

Again, the formula ( [5. 16| j easily follows from ( p.l| ). Also, a similar expression for T At (t) is 
obtained from ( |5.16|) by the replacement q — > g _1 . 

The proof of the basic relations (|5.1|) and ( j5.2|) given in Appendix is a straight- 
forward (but rather more technically complicated) generalisation of the corresponding 
Ug(sl(2)) results of JiJ. We consider more general T-operators which correspond to the 
infinite dimensional highest weight representations of U q (gl(3)). These new T-operators 
are defined by the same formulae as ( [4 . 2 1| ) 

T+(f) = Tr w + [e^ Ph L+(t)], L+(t) = <(*)[£], (5.18) 

except the trace is now taken over the infinite dimensional representations 7r+ of U q (gl(3)) 
with arbitrary highest weights /i = (/ii, /i 2 , ^3) the highest weight vector | 0) defined by 



( |2.41|) . Similar modifications to ([4.36|) allow us to define the operators T^ (t) related 
to the second evaluation map fl4.18|) . Then, by considering the tensor product of three 
(/-oscillator algebras corresponding to the product of three Q»(t) operators, we show that 

T /U,,aW = ^o 1 Qi(^ 2A11+2 )Q 2 (^ 2At2 )Q3(^ 3 - 2 ), (5.19) 

and similarly 

= -^ 1 Qi(^-^- 2 )Q 2 (^ 2 ' 12 )Q3(^ 2w+2 )- (5.20) 

Then, we make use of the Bernstein-Gel'fand-Gerfand (BGG) resolution of the finite- 
dimensional modules. The BGG result allows one to express the finite dimensional highest 
weight modules in terms of a direct sum of infinite dimensional highest weight modules. 
This implies that the T-operators for a finite dimensional module can be written in terms 
of |TTJ), ( pq ) as 

T,(t) = E(- 1 )' MT ^ + P)-pW' ( 5 - 21 ) 

<x6S 3 

where p = (1,0, —1), the summation is over all permutations o G S3 of three elements, 
and 1(a) is the parity of the permutation. There is an analogous result for the operators 
T Mli ^ 2) ^ 3 (t). The determinant formulas ( |5.1| ) and ( |5.2j ) follow immediately. 
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6 Conjectures: Exact Asymptotic Expansions 



6.1 Assumptions 

We will make the following assumptions 

i) Analyticity. The operators Aj(£) and A;(i) (and hence all the T-operators intro- 
duced above) are entire functions of the complex variable t in the sense that all their matrix 
elements and eigenvalues are entire functions of this variable. Their leading asymptotics 
at large t are given by 



logAi(t) ~ M(+it)* 



logAi(t) ~ M(-it) 



3(1-9) 



OO 



(6.1) 



where M is some c-number constant depending on the central charge c only. 

ii) Location of zeroes. For real p and large t the zeros of the eigenvalues of Ai(t) 
and Ai(t) accumulate along the positive imaginary and negative imaginary axes of t 
respectively. The asymptotics ( |6.1|) are valid everywhere except in the vicinity of the 
lines of zeroes (i.e. when | arg(it)| < 7r for Aj(i) and when | arg(— it)\ < 7r for Aj(t)). 

As is known the zeroes of the eigenvalues of Q-operators satisfy the Bethe Ansatz 
equations. In our case there are six sets of these equations involving different pairs of 
eigenvalues {Ai(t),Aj(t)} with % ^ j. To simplify the notation let us choose a particular 
pair of eigenvalues, say {Ai(t), A 3 (t)}, and denote by {t k }\^ =1 and their sets of 



zeroes 



Ax(t k ) = 0, A 3 (t k ) = 0, (6.2) 

ordered as \t%\ < \t 2 \ < ... and similarly t k . The equations (|5.10|) , ( p.llj ) (together with 
the standard arguments based on the analyticity of T-operators) imply the following 
Bethe- Ansatz equations 

'A^hq 2 ) Mhq- 1 ) 



7T- lo § 

2m 



2m 



log 



lMt k q-*) 
A 3 (t k q 2 ) 



A 3 (t k q) 
Atihq- 1 ) 



A 3 (t k q- 2 ) A x $ k q) 



-2e 3 p 



n 



(i) 



-2exp -4 3) + 1 



(6.3) 



2 ■ 



where n 



(1) „(3) 



G Z. Different eigenvalues of the operators Ai(t) and A 3 (t) correspond 



to different sets of integer phases {n 



(1) „(3) 



n, 



Jk=V 



The following assumption concerns the structure of the vacuum eigenvalues A^ ac \t) 



and A 3 vac \t) of operators Ai(t) and A 3 (t) acting in JF P . 

Hi). For real values of 2{wip) < g and 2(w 3 p) > —g all zeroes of A\(t) ( A 3 (t) ) 
located on the positive (negative) imaginary axis of t and the integer phases in (RO) run 
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over consecutive positive integers 

nk )=n k ) = k > fe = l,2,...,oo. (6.4) 



6.2 Results 

Following || consider the following operator valued function 

ll J r(i^/3)r(-zv(^ + i)/3)r((-i + zv)/3)r((i + zv)/3) y_ i00 t 1 J l[h - 

(6.5) 

where the integration contour goes along negative imaginary axis and 

? = T A_. (6.6) 

This integral converges for 3(g — l) < Sm v < — 1, however the definition of ^j(^) can be 
extended to the whole complex plane of v by analytic continuation. For example, using 
the product representation which holds for any eigenvalue Aj(t) of Aj(£) 

oo 

^W = n( 1 -sr). < 67 > 

n=l r ™ 

one can write the corresponding eigenvalue of as a Dirichlet series 

/ \ -if (1+0 

(r(i-^)j r(u/(f + 1)/3) » , M1+0/3 
^' (z/) ~ r(z^/3)r((-i + ^/3)r((i + z,)/3) ^' lt - ] w ' (6 ' 8) 

For < g < 2/3 this series converges absolutely for Qmu < — 1 and can be analytically 
continued to the whole complex //-plane by means of the standard technique f5"6|| . We 
expect that the function thus defined enjoys the following analytic properties 

Conjecture 1 The function is an entire function of the variable v. 

Converting the integral transform in (|6.5| ) one obtains 

i og A,(t) = -L/ -rK/3)r(i-Me + i)/3)r((-i + ^)/3)x 

2tt« Jc v v 

T((l + iv)/3) (T(l - g^"^ (it) iu{1+ ^ 3 9(y), (6.9) 
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where the integration contour C v goes along the line Sm v = — 1 — e with arbitrary small 
positive e. 

The values of at special points on the imaginary axis of v (where the Gamma- 

functions display poles) are of particular interest. For example, using (|6.8|) , (|6.7|) and the 
definition of the nonlocal IM Hn' in (|4.55|) it is easy to see that 

* 7 (-3m(l -g)) = —, — ; — K — T I — N , , 7 tt H^. (6.10) 

1 y)) T(ng)T(-l/3 + n(l-g))T(l/3 + n(l-g)) n 1 ; 

The other special values of of the imaginary axis are determined by the following 

Conjecture 2 The function has the following special values on the positive imagi- 

nary axis 

*,.((2n-l)z) = ^-( — ) 3 J g n Ian-i, n = 0,1,2,... (6.11) 

m = ^vrAf) vfr^ !rU - n=1 ' 2 '- (6 - 12) 

^ ain^" 1 -! ) = _ ,\ ; -mr/^ 71 > n=l,2,... 

r(n# x )r(-3+n(i-p 1 ))r(3+n(i-^ x )) 

(6.13) 

*.(0) = (6.14) 
J 2^^ 

where j = 1,2,3 and I_i = I is the identity operator. 



6V3tt V 3 J r 



2 

2 fx2n m - 



Below it will be more convenient to work with the variable s introduced in ( 4.56 ) 
instead of the spectral parameter t and exhibit all arguments of A(s) 

oo 

Afagtftu)} = exp ( - j^HgV 1 ). (6.15) 



71=1 



The above conjectures allow us to derive the asymptotic expansion of the operators 
Aj(t) for large t. 

Calculating the integral ( |6.9| ) as a formal sum over residues in the upper half-plane 
3m v > — 1 one has 

A,{s, g, <j>{u)} ~ C,{s, 0(it)} a"— exp { ^ 5 fc S "*/3(i- fl ) l fc J Aj { S ", g~\ g^u)}, 

k=— 1 

(6.16) 
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where 



oo 

Ajis-^g- 1 , g- l <t>{u)} = exp ( -^H« S -?), (6.17) 



n=l 

and the coefficients B n read 



= .3(1-^3(1-^ ^ n = 0,1,2,... (6.18) 

3(l-q)n\T 



= .3 ( i- g) ; V3(i- g -v ^ n = 12 (fU9) 

(l- 9 )n!r(--) 

An exact form for the operator Cj{g, <fi{u)} is unknown. However, an asymptotic vacuum 
eigenvalue in a special case can be calculated (see ( |C.56| )). Note, in particular, that the 
coefficient M in the leading asymptotics (|6.16|) is given by 

^feMs)^.^ (620) 

Naturally, considerations similar to those above are valid for the operators Aj(t) as 
well. The corresponding formulae are obtained from ( p.9|) simply by a replacement of 
Aj(t), tyj(v), Hn^ by their "barred" counterparts 

Aj(t) -> ~Aj(t), - Hj? - Hi j) , (6.21) 

accompanied by sign changes 

Wj -> -IB,-, I 2n -> -Ian, (6.22) 

where j = 1, 2, 3 and n — 1, 2, . . . , oo (the sign of the "odd" local IM l2n-i remains the 
same). Also, in the integrand of ( |6.9| ), the change t — > — t needs to be made. 

We can use these asymptotic results for the operators Aj(t), Aj(£) to calculate an 
asymptotic expansion for T(£). Using ( 5.1 0|) , we can write T(i) in the following form 

T(t) = Mqh) +Ar 1 (q-h)A j (qh) + AjV 1 *), i^j, (6.23) 

where 

Aj(t) = Aj = (6.24) 

jW Q,(^)' J QjCr 1 *) 
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Then, using ( |6.9|) we have 

log a, = 2„KP) + / ^r ( i-Mf + i)/3)r((-i + « ; )/3)r((i + fe )/3) x 

Jc v » r(i- «i/^/3) 

(r(l - ^)) W(1+ °(^) w(1+ « )/3 ^(i/), (6.25) 
which gives via the sums of residues in the upper half-plane v > — 1 

oo 

logAj(t) ~iml(jt)^ + J^C' fe (it) 5 ^J I fc , t -> +oo (6.26) 

fc=i 

where 

ttz = l3 , 3 9 J ( r(l - #) ) 1-9 , (6.27) 

r(i- |)r(|) v 1 y, J 1 ; 

and 

C k = -2ig k ^ sin (^) (r(l - g))'^ B k . (6.28) 
A similar expression is obtained for logAj(i), using (|6.21|) , (|6.22|) . It is easy then to see 



that only the second term in (|6.24j ) will contribute in the asymptotic limit, as the other 



two will be exponentially small. Thus we have 

oo 

logT(t) ~ mt^iTI-2^C 2n cos(^— J (t) W^)\ 2n 

n=l 

+ 2iJ2 C 2 „-i sin C K{2n ~ 1) ) (ty^hn-i. (6.29) 

n=l ^ ' 

The corresponding formula for T(t) is obtained from the above one by the replacement 
(|6~22|) . This result can be regarded as the the quantum version of (|2.29|) . In fact, if we 
take the limit £ — > (equivalently q — > 1), then, using the correspondence ( |2.20| ), we 
recover exactly the classical expression (|2.29|) .p| 



As in [[|] the main motivation to the conjectures given above came from the study of 
the large p asymptotics of the vacuum eigenvalues of the operators Aj(t) (see Appendix 
Q) as well as consideration of the classical limit g — > in the T-operator ( |6.23D . In this 
way we obtain 



2~( - — ^ 



¥f(f)|p-«x> = ; — 2-Fi( — + —,- + —,- + —, 1 + — ^)|,-th, (6.30) 

A )lp v / 3Tr(l/2 + W3) U 6 6'6 6'2 3' Agl )b 1 ; 



3 Recall the correspondence of variables t = A 3 
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where the subscript j—th means that an appropriate branch of the hypergeometric func- 

(6.31) 



tion is chosen such that (|6.14|) holds, and g 2 , g 3 are defined by 

g 2 = XiX 2 + X1X3 + X 2 X 3 , g 3 = -£1X2X3, 

where 



faP), J = 1,2,3. (6.32) 

The expression ( |6.30|) together with ( |6.9| ) allows one to derive the asymptotic expansion 
(|6.16|) in the limit p — > 00 which provided a basis for the conjectures ( |6.11|) , (|6.12|) , (|6.13| ), 

(EH)- 

The formulae ( |6.16|) allow us to obtain asymptotic expansions for the zeroes of the 
eigenvalues of the operators Aj(t) and Aj(t). Let us introduce new (rescaled) variables 
Eu\ Eu \ k = 1,2 (whose meaning will become clear in Section 7.2) to denote the 



J k ' k i 

positions of these zeroes 



where 



E k - % 9 t k J 



fc = l,2,..., 



Substituting ( |6.16| ) into the Bethe Ansatz equations one obtains 

00 



k k 2 g 



p» J « [ E - 



re=-l 



J k 



n=l 



OO re OO « .» 



k — n k 2 







n=-l 



n=l 



OO (6.33) 

(6.34) 

n 

% (6.35) 
ff , (6.36) 



where 



/3 n = -7T- 1 3 n <f "A 1 ^ sin(7rn/3(l - g)) B n , lk = tt" 1 3 3{1 ~ 9)n/g g 3n wn(im/g), 

(6.37) 

~ ~ (i) 

the coefficients £>„ are given in (|6.18|) , (|6.19|) and #4 and i? n are the corresponding 
eigenvalues of the local IM ( J2.8| ) and the dual nonlocal IM, defined in ( 4.7(J| ). Inverting 
the series ( |6.35| ), ( |6.36| ) one can express Ejp and E^ as asymptotic series 



E { k j) ~ T(AT ; (i) 



4 J) ^ T(ivi j) ), 



(6.38) 



where, for example, for g < 3/5 one has 



T (2) ~ a x 3 " 39 (l + a 2 x~ 2 + a 3 x~ 3 + 0(x~ 4 ) + 0(x 2_3/9 )) , x -> 00 (6.39) 
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T(x) ~ a x 



3-3 9 ^ + 



a 2 x~ 2 - a 3 x~ 3 + 0(x" 4 ) + 0(x 2 ~ 3/9 ) ) , 



X — ► oo 



«o = (^i) 3(3 - 1} 



"2 = il 



3T(-)r((e+l)/3) 

^r(|) 



3(l-</)/i 



-1 3(1-9) 



«3 = /?2/? 2 1 ^2 



7r(cot(vr(£ + 1) /3) - cot tt/3) ' 



3 5 3(1 



r 



2r ( ]) 3 r 



|(e + i))r((e + i)/3 N 



(6.40) 

(6.41) 
(6.42) 
(6.43) 



The Bethe Ansatz equations ( |6.3| ) can be solved numerically. We have performed 
such calculations for vacuum eigenvalues of Ai(t) and As(t) for a range of different values 
of g and p and compared the results with the asymptotic expansions (|6.38| ). A typical 
example of these calculations with 



1 



P 



1 _3_ 

2' 10 



(6.44) 



is presented in Table [l] where the numerical values for E's are compared with the corre- 
sponding values given by six terms of the asymptotic series ( |6.38| ) (the value of the leading 
term in ( |6.38|) is also shown). 



6.3 Special cases 

The following special cases are related to a particular choice of the vacuum parameter 
p = (p, 0) such that 



(w 1 p)=p, (w 2 p)=0, (ttf 3 p) 
and the weights ( |3.24|) read 



-P, 



A 2 = — + 



2 c-2 



A, = 0. 



(6.45) 



(6.46) 



9 24 

The vacuum eigenvalues in this case obey the symmetry 

A[ mc \t) = A { r\~t), 4 mc) (t) = 4 mc) H), 4 mc) (*) = A { r\-t), (6.47) 



T {wic) (t) = T (vac \-t), 



(6.48) 
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n 


4 1] (BA) 


^i 1} (AS) 


£i 1} (LT) 


Ei 3) (BA) 


(AS) 


Ef(LT) 


1 


426.44998 


429.22624 


615.63 


42.72214 


25.28503 


147.89 


2 


871.85882 


872.59835 


1069.31 


232.30688 


229.38814 


371.98 


3 


1466.71208 


1466.98579 


1674.60 


563.02921 


562.24812 


723.05 


4 


2224.92354 


2225.04588 


2443.32 


1039.72260 


1039.43412 


1216.12 


5 


3157.57647 


3157.63845 


3386.10 


1673.43618 


1673.30798 


1864.09 


6 


4274.50955 


4274.54393 


4512.70 


2474.97610 


2474.91149 


2678.44 


7 


5584.74586 


5584.76628 


5832.14 


3454.37424 


3454.33857 


3669.55 


8 


7096.68686 


7096.69967 


7352.86 


4620.94608 


4620.92497 


4846.97 


9 


8818.22784 


8818.23622 


9082.79 


5983.39245 


5983.37925 


6219.54 


10 


10756.83708 


10756.84278 


11029.43 


7549.88941 


7549.88079 


7795.56 


11 


12919.61477 


12919.61876 


13199.92 


9328.16145 


9328.15561 


9582.83 


12 


15313.33879 


15313.34165 


15601.06 


11325.54050 


11325.53642 


11588.77 


13 


17944.50152 


17944.50363 


18239.38 


13549.01388 


13549.01096 


13820.41 


14 


20819.34016 


20819.34174 


21121.13 


16005.26365 


16005.26150 


16284.47 


15 


23943.86187 


23943.86307 


24252.34 


18700.69926 


18700.69765 


18987.41 


16 


27323.86522 


27323.86615 


27638.82 


21641.48502 


21641.48379 


21935.42 


17 


30964.95845 


30964.95918 


31286.21 


24833.56334 


24833.56239 


25134.47 


18 


34872.57528 


34872.57586 


35199.94 


28282.67468 


28282.67393 


28590.32 


19 


39051.98867 


39051.98914 


39385.30 


31994.37469 


31994.37410 


32308.54 


20 


43508.32291 


43508.32328 


43847.43 


35974.04923 


35974.04875 


36294.55 


21 


48246.56426 


48246.56456 


48591.33 


40226.92739 


40226.92700 


40553.57 


22 


53271.57046 


53271.57072 


53621.85 


44758.09307 


44758.09276 


45090.72 


23 


58588.07919 


58588.07940 


58943.76 


49572.49523 


49572.49497 


49910.94 


24 


64200.71562 


64200.71580 


64561.67 


54674.95691 


54674.95669 


55019.08 


25 


70113.99930 


70113.99945 


70480.11 


60070.18348 


60070.18330 


60419.85 



Table 1: 
flOP for 

4° (AS) 



Comparison of the numerical solution of the Bethe Ansatz Equations 
En \ En^ with g, p given by (|6.44| ), to the value from asymptotic series (|6 
and the leading term of the series En (LT). 
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while the expression ( |4.81|) reduces to 



H (i)(vac) = ^ 3 r(l-f) T(g + p)T(g/2 + p) 

2sm(7rg)T(l-l)T(l-g)T(l-g + p)T(l-g/2 + P y l ' ' 



Further the asymptotic expansion of the vacuum eigenvalues T^ ac ^ (t) given by fl6.29| ) now 
reads 

oo _ 

T^)(t) =mtW^ + 2^cos((n-2)7r/3)C' 2 „_ 1 /£ a _ c it^, | arg(t) | < \^~9)- 

n=l 

(6.50) 

6.3.1 The case g = 1/2, c = -10 
Consider the case when 

c=-10, 9 = 1/2, q = e i7T9 = i. (6.51) 



Using ( |4.57|) and (|6.47|) it easy to see that the Bethe Ansatz equations ( |6.3| ) in this 



case reduce to those for the t/ g (s/ 2 ) related Q-operators studied in ||. This fact has 
been noticed in [^(J and recently discussed in |57]]. Therefore the functions A^ ac \t) and 
A^ ac \t) coincide (up to a normalisation of the spectral parameter t) with A+ ac \t) and 
y^vac) ^ f rom re f || |7j T/he re l eV ant parameters used therein are related to our notation 
as follows 

c* = ~ P 2 = \, V sh = % A = ^= P 2 -9/16. (6.52) 



Comparing now the asymptotic expansion ( |6.16|) with the corresponding result of ref 



(the eq.(4.19) therein) and using the expressions ( [4.73D and ( |6.49|) one can reproduce the 
values of first three local IM j^ vac \A, = -25/2) given in eq.(40) of ref . @ 

/(^)(~ C )(A )C ^) = l/f ac) (A 2 ,0,c) = A + H|, (6.53) 

2 48 

4 s ' 2)(mc) (A,c^) = -^H[ 1 ^ ac \l/g,p/g)^A 2 + 7 -A+^, (6.54) 

/| S,2)W (A ; ,c sh ) = -li vac) (A 2 ,0,c) = A 3 + -A 2 - ^^A + (6.55) 
s v ' ; 8 5 V 7 16 98304 774144 v ; 



7 This correspondence extends to certain excited states as well (see also Section 8.2), but we will not 
elaborate this point here 
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6.3.2 The case g = 2/3, p = 1/6 

Consider now the vacuum eigenvalue ( |6.50| ) of the operator T(t) when 

c=-2, £ = 2/3, p = 1/6. (6.56) 

The value g = 2/3 does lie in the domain Q3.38] ) therefore the results of the previous 
Sections do not directly apply. In particular the definitions ( |4.31| )- fl4\34j) require a renor- 
malization since the nonlocal IM ( [4.33j ) diverge logarithmically at g = 2/3. To obtain the 
asymptotic of T(t) at large t in this case we set g = 2/3 — e in ( |6.50| ) and let e — > (this 
corresponds to an analytic regularization of divergent integrals). The first few terms of 
the formula (|6.50|) then read 

l ogT (- C ) (t) | g=2Ap=1/6 = -V3atlogt + Ct+ ^(at)- 1 - ^p(^)~ 5 + • • • (6-57) 

where C is a renormalization constant and 

a = T(l -g) 3 = T(l/3) 3 . (6.58) 

On the other hand the eigenvalues of T(t) can be computed independently. Using the 
results of one can show that in the case ( |6.56| ) the fusion equations ( |5.14| ) for the 
eigenvalues of T(t) simplify to a single equation 

T(tq*)T(tq~*) = 2cosh(a7rt) T(t), q = e 2m/3 (6.59) 

which with an assumption that T(t) is an entire function of t and has the leading asymp- 
totics at large t as in ( |6.57| ) can be completely solved 

m , \ ™,f, M rW n A (1 + 2n fc - 2ate^/ 6 )(l + 2n k - 2ate"^/ 6 ) 

T(t) = T (mc) (t) TT ) -j^, 6.60 

w Wil l + 2n h + 2ate™/ 6 1 + 2m + 2ate"^/ 6 V ; 



fc=i 



{> r(i + ate-/6)r(i + ate— /6)- ^ Di ^ 

where the integer L > 0; {n^} is an arbitrary increasing sequence of integers < rii < 
n 2 < ■ ■ ■ < Hi an d C is an arbitrary constant. The vacuum eig envalue T( mc )(t) is a 
unique solution of ( p.59| ) with no zeroes in the strip |arg(t)| < n/3. One can easily check 
that its asymptotic expansion at large t perfectly agrees with ( |6.57|) . 



44 



6.3.3 The case g = 5/6 

Below we will compare the asymptotic expansion ( 6.50|) for g = 5/6 against the numerical 



results of [58[ for two values of p 



p=l/6, p = l/3, 



(6.62) 



which corresponds to 

A 2 = 0, A 2 = 1/10, (6.63) 

respectively. Again, for g = 5/6 the formula ( |6.50|) requires a renormalization. Using the 
analytic regularization as in the previous subsection one obtains 



T {vac \t)\ 9=5/6 = ~{atf log* + Ct - ^5 (at)- 2 l[ vac) + 3 ™lf (at)- 10 J 5 (wc) + . . . 
v3 



43472 



where 



a = r(i-s) 3 = r(i/6) J 



(6.64) 
(6.65) 



and I^ ac > and J^ vac) are given in ( |4.73| ). 



r(vac) 



ijTHc, A 2 ,A 3 ) 


Numerical values |p8 


if ac) (6/5,0,0) =-i =-0.05 


-0.499999999599 10" 1 


lf ac) (6/5,0,0) - ^ - 0.115238095238... 10" 1 


-0.115238095240 10" 1 


l[ vac) (6/5, 1/10,0) = ^ =0.05 


0.500000000002 10" 1 


lf ac) (6/5, 1/10,0) - 4 2oo - 0.497619047619 ... 10~ 2 


0.497619047634 10~ 2 



Table 2: Exact values of I^ vac \c, A 2 , A 3 ) for g = 5/6, and p given by (|6.62| ), compared 
with the corresponding numerical results from 



As shown in the eigenvalues of T(t) in this case satisfy the equation 

T(tg)T(tg- 1 ) = e^r' 2 T(-t) + e"^ 2 T(t). (6.66) 

In this equation was transformed to TBA-like integral equations, due to Kliimper and 
Pearce ||59|| , and studied numerically for a few vacuum eigenvalues of T(i). In particular, 
the numerical values of the coefficients in (]6.64|) were calculated. We have verified that 
the values of the local IM extracted from this numerical data are in a good agreement 
with the corresponding exact values determined by (|4.73|) (see Table pj). 
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7 Applications 



7.1 Boundary Affine Toda Theory 

The Hamiltonian corresponding to the (non-equilibrium) boundary affine Toda theory 
(BAT) with zero bulk mass reads 

Hbat = Hq + Hi, (7.1) 

Ho = ^ f dx(U 2 + ), H x = - A g Ut e i(e,,#B+v<) j (7 2) 

where g>, k and V = (Vi, V 2 ) are parameters; <fr(x) = ($i(x), $2(0;)), U(x) = (IIi(x), U 2 (x)) 
are field operators obeying canonical commutation relations 

[U a (x) , ®b(x')} = -27rig 8(x - x')S ab , a, 6 = 1,2 (7.3) 

the variables ui, u 2 , u 3 

u 1 u 2 = qu 2 u 1 , u 2 u 3 = qu 3 u 2 , u ?J Ui = quiU 3 , u x u 2 u 3 = , (7.4) 



describe boundary degree of freedom and = $(0). The way g enters the Hamiltonian 
(771) allows one to interpret it as a quantum parameter, because it always appears in the 



combination gh; in what follows we will set h = 1. We will also assume that the parameter 
V satisfies the requirement 

( ei V)>0, (e 2 F)>0. (7.5) 



At a nonzero V and a temperature T the system ( |7.1|) develops a stationary non- 
equilibrium state which can be thought of as the result of an infinite time evolution of the 
equilibrium state of the corresponding "free" system, with the interaction term (the last 
term in ( |7. 1| ) ) adiabatically switched on. We will denote by ( A)n the expectation value 
of an observable A over this non-equilibrium stationary state. 

The density matrix P(t) of the system is determined by the Schrodinger equation 

d t P(t) = -i[U(t),P], (7.6) 

with the initial condition 

P(t)| t= _oo = P = Z, 1 e~ RHo , R = g/T, (7.7) 
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where Po is the equilibrium density matrix of the free system at the temperature T. Using 
the interaction representation corresponding to the Hamiltonian H in (|7.1|) one can write 
the density matrix P(t) as 



P(t) = e - iHot S(t, -oo) P S(-oo, t) e iHot , (7.8) 



where 



S(Mo) = rexpj-^y dr H? nt) (<r)} = 

oo „ t 

= 1 + / %({r})Hf n * ) (r 1 )H? nt) (7i 1 )...Hi i ^(T B ) ) 

./in 



(7.9) 



with 



H ( l nt \t) = e iHot Hi e~ mot . (7.10) 
In ( [7.91) and below the shorthand notation for the multiple ordered integrals 

/ Vk({r})= dn dr 2 --- dr k . (7.11) 

J to J to J t() J to 

is used. The expectation value of an arbitrary operator A over the stationary non- 
equilibrium state ( |7.8| ) (which will be marked with the subscript N) reads 

(A) N = Tr H [P(t)A]. (7.12) 

The trace is taken over the space 7i = T ®U where T is the Fock space representing the 
Bose commutation relations ( |7.3| ) and U is the space of states of the boundary variables 
Ui, % = 1,2,3 (with the trace over these states defined exactly as in ( |4.60| )). Substituting 
(O) into ( [77T2D one gets 

(A) N = Tr w [P S(-oo, t) A^(t) S(t, -oo) ] = ( S(-oo, t) A^ mt \t) S(t, -oo) ) , (7.13) 

where (. . .) denotes the expectation value over the equilibrium state ( |7.7| ) of the free 
system, and the superscript "(mt)" means that this operator is taken in the interaction 
representation, i.e. PsS mt \t) = e !H(|i Ae _!Hl)1 . Equivalently, one may write the above 
expectation value as 

<A) N = Tr w [PA(t)] , (7.14) 

where P stands for density matrix of the system at t — 0, i.e. P = P(0), and A(i) is the 
full Heisenberg operator 

A(t) = S(0, t) A (mt) (t) S(t, 0) . (7.15) 
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All above formulae are very well known (see e.g. fl6"0|); the notation used here is exactly 
the same as in |HJ. We are interested in the expectation values (Yj )n of the Heisenberg 
operators 

Yj(t) = S(0,t) Uj vf nt \t)S(t,0) . (7.16) 

where 

Vf nt \t) = e *y.*2 nt) (*)+vt) , ■ = ^ 2? 3 {7 U) 

and &g nt \t) is the boundary field &b in the interaction representation. It is convenient 
to introduce also auxiliary operators 

V,(Mo) = S(t ,t) Uj Vf nt \t)S(t,t ) , (7.18) 

where to is a parameter. For to = ([7.18Q coincides with the Heisenberg operators (|7.16j ), 
and according to ( |7.13| ) the expectation values of ( |7.16| ) can be expressed through ( |7.18| ) 
as follows 

(V ff (t)) N = lim (Y a {t,t )) , (7.19) 

to— >— oo 

where a = 1, 2, 3. 

Using the series expansion for the evolution operators ( |7.9| ) and the commutation 
relations 

Vf nt \h)Vf nt \t 2 ) = q 2 ^Y { ; nt) {t 2 )Yt nt \ti), h > t 2 (7.20) 
one can represent (|7.18| ) as a series of time-ordered integrals of products of \^ nt ^ 

Y a {t,to)=uM nt \t) 

+V (int) {t) J2 C k (a,a u ...,a k ) / V k ({t}) V^\ tl ) ■ ■ ■V^ t \t k ) , 

k = l <7l,...,<7 fe = l,2,3 Jt ° 

(7.21) 

where the sum is taken over all arrangements of the "charges" <7i, . . . , a k = 1, 2, 3. The 



coefficients C k (a\ <j\, . . . , a k ) can be easily calculated following the arguments of ref. |30 
From the definition ( |7.18| ) it follows that 

i J-V,(Mo) = [HS m4) (t ),V ff (t,t )]. (7.22) 



Substituting the the series ( |7.21| ) into this differential equation one gets a recurrence 
relation for the coefficients Cfc(cx| a\, . . . , <J k ), which (with the initial condition Co(c) = u a ) 
has a unique solution 

C k (a, a 1 ,...,a k )=u i r Xk W WlXk) (—)" q F{Xk) f[sm(ng w^X^) , (7.23) 

9 3=1 
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where 

j 

(jj = a, j - 1 (mod 3), Xj — ^ e as , a = a (7.24) 

s=0 

and 

1 / \ 2 1 / \ 2 1 

FX): 1 " N 



2 



(ti^x) -2^) (7-25) 



Using ( |7.21]) , (]7.23|) one obtains an infinite series representation for the expectation 
values (|7TT9D F1 

<V <T ) N = : V Ml ^^...^Sn-lP , 

k sin ire ^ — f z — ' V - LJ - sin iro / 

n=l <Tl,...,CT3n-l i = i v ^ 

(7.26) 

J(a,a 1 ,...,a^ 1 \p)= f V^r}^^ W ]J (2 sin (Z^p)) 9 ^. 

(7.27) 

where cr, = 1,2,3 and the internal sum taken over all possible of values of {0^}, i = 
1, . . . , 3n — 1 such that 

e ao + e ai + . . . + = . (7.28) 

It is convenient to represent the sequence of the "charges" 00 > °i ; • • • > a k f° r each term in 
(|7.21|) by a path staring from the origin X = of the triangular lattice 



X = met + n 2 e 2 , ni,n 2 eZ, (7.29) 
and formed by the sequence of the "step" vector e CT0 , e U2 , . . . , e CTfc . The vectors Xj defined 



in (|7.24| ) then represent the end point of the path after j + 1 steps. It is important to 
note, that not every such path contributes to ( |7.21| ). For example, for <7o = 1 only paths 
which are 

i) always contained in the sector (see Fig [j]) 

WiX m > 0, w 2 X m >0, m = 0,l,...k (7.30) 

ii) do not return to the origin (except for the final step), i.e., 

X m ^0, m = 0, 1 (7.31) 



Representations of this type for a similar problem related to the boundary sine-Gordon model were 
first derived in H] by combinatorial methods. 
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correspond to nonzero coefficients Ck(ao, a±, . . . , o^), since otherwise one of the factors 
smijig w a i.Xj_i) in ( |7.23| ) vanishes. Obviously, the paths (|7.28|) contributing to ( |7.26|) 
return to the origin at the last step. As a result all integrals appearing in ( |7.26|) converge 
at large r's for 

2(eip) > -g, 2(e 2 p) < g. (7.32) 




Figure 1: The shaded area is the region defined by ( |7.30|) . We also show an admissible 
path X 5 = ei + e 1 + e 3 + e 3 + e 2 + e 2 . 



The expectation values (|7.26 ) coincide with equilibrium expectation values of certain 
operators in the system which is similar to Q7.1J ) but contains different boundary degrees 
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of freedom obeying the g-oscillator algebra ( 4.43| ) . Consider the Hamiltonian 



H = 4^ dx ( n2 + *l ) + \ hv ~ y 9 E * eie '* B ■ ( 7 - 33 ) 

where «fr(x), U(x) are again the Bose field operators obeying the same commutation 
relations ( [7.3| ) 

2 

ft = -(/11C1 + h 2 e 2 + h 3 e 3 ), (7.34) 

and 

V\ = (q-q" 1 )^, V2 = (q- q^afq* 2 , y 3 = q?(q - g -1 ) 2 g- Wi, ai a£, 
/ii = — 7ii — 27i 2 , ft<2 = 27i 1 +7i 2 , /ix = —7^1+7^2, 

(7.35) 

where the operators a^ and 7ij, i = 1,2, commute with 3>(x), U(x) and satisfy the 
commutation relations Q4.43 ). Let p be any representation of the algebra ( |4.43|) such that 



the spectra of Tii and 7i 2 are both real and bounded from above. The Hamiltonian ( fT.l|) 
acts in the space 

H = T®p, (7.36) 



where T is the Fock space representing the Bose commutation relations (|7.3| ) and for the 
values of V satisfying (|7.5|) this Hamiltonian is bounded from below. Then the system 
( |7. 1| ) has a thermal equilibrium state described by the standard density matrix 

P = Z(k, V)- x e- m , R = g/T (7.37) 

where 

Z(k, V) = Tr^e"™, (7.38) 

is the partition function. Denote by (A)e the expectation value of an observable A over 
this thermal equilibrium state 

(A) E = Tr^(AP). (7.39) 
We will be particularly interested in the expectation value (Vi)^ of the operator 

V 1 = y 1 e i «*», (7.40) 

and will show below that 

(Vi> B = (Vi>jv, (7.41) 
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where (Vi)jv is given by (|7.26|) above. 

The proof is closely parallel to that used in the proof of the analogous result for the 



Boundary Sine-Gordon model, found in [3C]. We will use the interaction representation 
with the Hamiltonian 



— 1 
H 



Aug 



° dx (n 2 + *l) + ±hV, 



(7.42) 



corresponding to k = in ( 7.331 ). Consider an auxiliary operator 

Vi(t, to) = S(t , t) e m ° l ?! e- iHo *S(t, to), (7.43) 
where S(t,t ) is defined as in ( |7.9|) but with the interaction Hamiltonian replaced by 



^9 i= i 



(7.44) 



(7.45) 



The operator ( 7.43|) can be represented in a form similar to (|7.21| ) 

V 1 (t,t ) = aZ Vf 4) (t) 

OO „£ 

+V (int) it) J2 C k (a,a u ...,a k ) / V k {{t}) V^(h) ■■ -Vg^fe 

fc = l <Tl,...,<7fc=l,2,3 

where 

C*(a, n,. ..,<*) = (aJ)^*)(aTg-" a )-^ (^)* [[M^ tn^-i) , 

(7.46) 

and the rest of notation is the same as in ( [7.21 ). In particular, the same remarks about 

the paths contributing to ( |7.21| ) apply to ( |7.45[ ). We can now use the representation ( |7.45| ) 

in (|7.39|) . Repeating the arguments which used in derivation of eq.(3.15) in ref. |3(J one 

obtains _ _ 

Tr H [ e -^Vi(0,t )] 



(V 1 ) B = lim 

Using now (|7.45|) and the simple property of the trace 

Tr^[e-^(ar) n (a 2 ) m ] = ^o^Tr^e^] 



(7.47) 



(7.48) 



one arrives exactly to the series ( [7.26 ) for the expectation value ( 7.471) . This proves the 
relation ( |7.41| ). 
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Further, comparing the Matsubara representation for the partition function ( [7.38 ) 



with the vacuum eigenvalues of the operators A*(t) defined by Q4.49Q one can conclude 

that 

5"' V) = A { ™ c \t,p), (7.49) 

z(o, v) 3 y ! y 1 

where the parameters t and p related to k, V and the temperature T as 

k sin Tig ( g \ -9 1 3 



t = i 



2ttT \2txT 



p = -i-i-V. (7.50) 



Expressing now the equilibrium expectation value (Vi)e through the partition function 



and using ( |7.41|) we obtain 



< Vx > E = ^TdJogZ(n, V), (7.51) 



(¥ 1 ) N = ^TdJogAl' ac \t,p). (7.52) 



In a similar way one can show that 



(V 2 ),v = 2 -Td K \ogA { r\t,p), (V 3 ) N = ^Td K logA^ ac \t,p). (7.53) 



Finally note the g — > 1/g duality properties of the expectation values ( |7.52| ), ( |7.53| ). 
Introduce the following combinations expectation values 

J a = inn ( V cr _i)jv — i^n (V^n, a = 1,2, 3 (7-54) 

then 

J a = 2mTtd t log f Y » ( 7 ' 55 ) 

where (a, ft, 7) is a cyclic permutation of (1,2,3). 

According to ( |6.16| ) the quantities ( |7.54j ) satisfy the following "strong weak" duality 



relations 

J a (g, V, -k) = -e„ V - g- l J a {g-\g V, -C(g)(-K)^), (7.56) 
C(g) = T -^( T -^)K (7.57) 



7T \ TC 
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7.2 Relation to third order differential equations 

Consider the differential equation 

{* + «i±M±M^ _ « + X 3«w (x) = (7 . 58) 

where 

£i+£ 2 + £ 3 = 3. (7.59) 

For 

3fte(A + 3) > Ue(£ 2 ), Ue{£ x + 3) > 3£e(4), (7.60) 
the equation (|7.58|) has a unique solution satisfying the condition 

ip(x, e, e u e 2 , e 3 ) = x £i + o(^ i+3 ), (7.61) 



This solution can be analytically continued outside the domain (7.60). Obviously the func- 



tions i/)(x, E, £ 2 , £1, £3) and if)(x, E, £ 3 , £\, £ 2 ) defined in this way satisfy the same equation 
(|7.58|) and for generic £i,£ 2 ,£ 3 the solutions 

*l> 1 (x)=i>(x,Ej 1 ,e 2 j 3 ), fo(x) =i>(x,Ej 2 ,e 1 j 3 ), fo(x) =i>(x,E,e 3 ,e l ,e 2 ), 

(7.62) 

are linearly independent since 

Wr [^{x)^ 2 {x) : ijj 3 {x)) = (h - £ 2 )(£ 2 - £ 3 ){£ 3 - £,), (7.63) 

where Wr denotes the usual Wronskian. 

Further, for all E the equation ( 7.58Q has a unique solution x{ x i E) which decays at 



x — > +00. We normalise this solution as 

x ( x ) = x -« exp | _ ^_ + 0(x a " 2 )}. (7.64) 

It can be expanded in the basis ( |7.62| ) 

x{x,E) = R x {E)^ x + R 2 {E)^ 2 + R 3 {E)^ 3 . (7.65) 

It was conjectured in |H| that the connection coefficients Ri(E) are simply related to the 
vacuum eigenvalues of the Q-operators (|4.49|). The exact correspondence is as follows: 



Ri(E) a (vac) 



Ri(o) 



Ar c> (t,p), i = l,2,3 (7.66) 
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where 



a 



1 

9 



E = i pt, 



3(pwi 



+ 1, 



(7.67) 



with the parameter p and the vectors Wi defined in ( 2.34f) and fl6.34|) respectively. The 
conjecture (|7.66|) was supported in |2(| by extensive numerical work, but no analytical 
proof is currently known (see Section 8.3 below for a discussion). 



E n 1] @ E n l) (AS) 



(1) 



2.88682816 
5.15225193 
7.13037330 
8.94036211 
10.6356087 
12.2451641 
13.7870634 
15.2734900 
16.7131734 



2.8787210 
5.1513633 
7.1300266 
8.9401849 
10.6355029 
12.2450946 
13.7870146 
15.2734540 
16.7131460 



p(3) 



PI E n 3) (AS) 



0.890944821 
3.643901890 
5.817136367 
7.737709470 
9.508485338 
11.17453999 
12.76111150 
14.28420813 
15.75481968 



0.66892875 
3.64527385 
5.81718647 
7.73772076 
9.50848916 
11.17454166 
12.76111235 
14.28420862 
15.75481998 



Table 3: Numerical values for E n an d Effl from Table 2 in ]2(|, truncated to 8 decimal places, 
compared with E^ (AS) and -E^^AS), obtained from our asymptotic formula, ( [7.71 ). 



Note that, due to the correspondence ( |7.66| ), the formulae ( |6.16| ) provide asymptotic 
expansions for the zeroes of the coefficients Ri(E). As an example consider the following 
particular case of (|7.58|) , 

d 3 x ^{x) + x = E^f(x), (7.68) 



which corresponds to 



and 



a = 1/3, 



g = 3/4, (pwx) = 1/4, (pw 2 ) = 0, (pw 3 



-1/4. 



(7.69) 



(7.70) 



The zeroes Ri(E) and Rs(E) in this case have been extensively studied numerically both 
from the solutions of the differential equation as well as from the non-linear integral equa- 
tion |18|-f20| which is equivalent to the Bethe Ansatz equations (|6.3|), (|6.4|) for the vacuum 
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eigenvalues A^^it). We have compared this numerical data against the asymptotic ex- 
pansion ( J6.38[ ). For g — 3/4 the expressions ( |6.35| ), ( |6.36j ) allow us to derive only one 
non-trivial term in ( |6.39|) , ( |6.40|) (note that the first dual nonlocal IM also contribute). 
For the values of p given in ( [7.70| ) the corresponding formulae read 



i-in. 



E, 



(3) 



/ 87T \3 
V373V 

/ 8tt \| 



in 



in 



_ 5>3/4 



1 + 256tt 2 ( n \} 2 
1- 15 



2567r 2 V"- §) 



+ 0{n~ A ) 



(7.71) 



The numerical values are presented in Table [3]. 

As shown in |52] the correspondence (|7.66|) allows us to derive a certain duality relation 
for the vacuum eigenvalues A\ vac \t,p). It is convenient to make a change of variables 
^(x) = e y ty(y), x = e y which brings the equation (|7.58|) to the form 



{dl - [v x v 2 + Vl u 3 + v 2 v 3 )d y + iv x v 2 v z + e 3 ^ y - Ee 3y ^(y) = 0, 



where 



Vj = —i(£j — 1) 



It is easy to see that the change of variables 

y = gy' + yo, 



-3yo 



-9 3 E, 



(7.72) 
(7.73) 

(7.74) 



transform ( |7.72j ) to an equation of the same form but with the parameters g, and 
replaced by 

g^g-\ v l =gv u E -> -^"f (-£)-?. (7.75) 

Using this fact in ( [7.651 ) one obtains the following duality relation for the connection 
coefficients FL 



e^# a (g,t/i,z/ 2 ,z/ 3 ,E) = R a (g \ gv Xl gv 2l gv 3 , -c(g) j-E) Q 

R P (g, V X , U 2 , U 3 , E) R^g-\g Vu g U2 , gu 3 , - C ( 9 ) (-E)' 



1 ? 



c(g) = g 3 ^, 



(7.76) 

where a, [3 = 1, 2, 3. Note that due to ( [7.66Q , this relation is in agreement with ( |6.16| ) and 

arm . 



8 Discussion 

In conclusion we make a few remarks on some open problems relevant to the topics of 
this paper. 
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8.1 Truncation of the functional relations 



It is well known that the fusion equation for the transfer matrices in solvable lattice models 
truncate when q 2 is a root unity. A similar truncation, of course takes place for the CFT 
T-operators. For the theories related to the quantum affine U q (sl(2)) this truncation is 
well understood. For example, consider the CFT T-operators Tj(A), j — 0, |, 1, | . . ., of 
ref. 0, associated with 2j + 1-dimensional representations of U q (sl(2)). If q = e lw/32 is 
a primitive root of 1 or —1, such that q N = ±1 with some integer N > 2 then one can 
easily show (see, e.g., H) that 

Tn(\) = 2cos(27riV» + Tn_ 1 (\) , (8.77) 

where p and j3 are related to the Virasoro highest weight A and the central charge c as 

A = |U^p c = 1 - 6(/3 - /T 1 ) 2 , ? = e^ 2 (8.78) 

As is shown in || this allows one to bring the fusion relations for Tj(A) to the form 
identical to the functional TBA equations (the ^-system) of D N type p9| , |6"3|[ . 
Further reduction occurs at particular values of the vacuum parameter p 

Tjv_i(A) = 0, when z N = ±1, z = e 2nip (8.79) 



which corresponds to the RSOS reduction in the six-vertex lattice model [64]. This latter 
RSOS-type reduction is also well known for the case of U q (sl(3)) Hl6| . For example, it 
readily follows from the determinant formulae ( |5.1| ), ( |5.2j ): when the vacuum parameters 
Zi defined in ( |3.31| ) become roots of unity zf = 1, i = 1, 2, 3 along with q N = 1 and either 
fii + 2 = iV or fi2 + 1 = iV the determinants (|5.1|) , ( |5.2|) vanish. However, no simple 
analogues of the reduction ( ^.77| ) valid for arbitrary values of Z{ is known (except for the 
particular case of Sect |6.3.1| which accidentally reduces to U q (sl(2))). Similar remarks 
apply to related problems such as the string structure of the Bethe Ansatz equations and 
the TBA equations. In particular, no U q (sl(3)) analogues of the type TBA equations 
are known. 

8.2 Connection with U q {A^ ] ) related IQFT 

The set of classical IM ( |1.2|) admits consistent reductions V{u) = or V{u) = U'{u)/2 
causing all even IM If, , k = 2,4 (mod 6) to vanish. Remarkably, the non- vanishing 
IM 4 d) , k = 1,5 (mod 6) remain in involution |T(J with respect to Virasoro Poisson 
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bracket of the field U(u). The quantization of these conserved quantities leads to new 
sets (with respect to the quantum KdV case @, ^3], ) of local IM in c < 1 CFT. 
Actually, there are two new non-equivalent sets of IM related to the different reductions 
mentioned above. This structure of the local IM naturally arises in the conformal limit 
of the massive IQFT associated with the twisted quantum affine U q (A^) |65|,|66i. The 
universal CFT L-operator for this case reads 

r r 27T 



e inPh Vexp 



e~ mu) : 2/o+ : e^ (w) : yAdu 



(8.80) 



where <fi{u) is the chiral Bose field normalised exactly as in and h, yo and y\ are the 



generators of the Borel subalgebra of U q (A^) where q is related to the Virasoro central 



charge c as in QS.78j ). In ref. |)7| Fioravanti, Ravanini and Stanishkov constructed an 
infinite set of finite-dimensional representations TT s (t), s = 0, 1,2, . . ., of this subalgebra 
(with the dimensions dim[7r s (t)] = (s + + 2)/2) and defined a set of T-operators 



T.(t) 



Tr Mt) [^ Ph C], 



S = 0,l,2, 



Here we use the spectral parameter t which is related to A of |)7J] as t = A 3 . These 
operators act in the Virasoro highest weight module Va with A defined in ( 3.78|) ; they 
are entire functions of t and their leading asymptotics at large t can be easily estimated 



logT s (t)~m s t 2 /( 6 - 3 ' 32 ), |t|-oo 



(8.82) 



where m s are some numerical constants. Further, at small t one has 

T x (t) = 1 + 2 cos(2vrp) + d t + 0(t 2 ), 
where the vacuum eigenvalue Gi in Va can be explicitly calculated 



(8.83) 



G 1 L( 2 ) = 2ir e 



2ir rui 

dui i du2- 

Jo 



e ip(ui+ " 2) (2 sin( 



2 sin — 2 sin — 
2 J \ 2 J 



47r 3 r(i + /5 2 /2)r(i - 3/3 2 /4)r(i - p 2 ) 



T(l + (3 2 /4)T{l + p- P 2 /2)T{1 -p- (3 2 /2)T{l +p- /3 2 /A)T{l -p- (3 2 /4) ' 
In |J7| it was shownQ that these operators satisfy the functional relation [^] 



T s (t g V2) Ts ( tg -V2) = Ts (_ t ) + T^T^it), q = e 



iwf3 2 



15) 



9 Although the calculations in |67| have been performed in the lowest non-trivial order in t the relation 
,'85| ) is expected to hold in all orders in t. 
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Now turn back to the W3 theory. Consider the W3 module Va 2 ,o, with the weight 
A 3 = 0. Let V sym be a subspace V sym G Va 2 ,o of this module invariant with respect to 
the action of the automorphism 013 defined in Sect. [4.1| Then from the definitions ( 4.21| ). 



(|4.36|) and the relations ( |4.18|) , (|4.11|) one concludes that the operators T m (t) and T m (t) 



from (|08D , (HD obey the relation 

T m (t) = T m {-t) (8.86) 

when acting in the subspace V sym . It follows then from (|6.29| ) and the discussion there that 
all even local IM k = 2,4 (mod 6), vanish in this subspace. Denote by T m (t\g,p)) 
the eigenvalues T m (t) in V sym G Va 2 ,o where g is related to the central charge c of the 
W3 algebra (|3.22j ) and p parameterises the weight A 2 = p 2 /g + (c — 2)/24. Taking into 



account flS.86|) the functional equations (|5.14| ) become equivalent to (|8.85| ) provided one 



makes a replacement q — > 52. Thus we have good reasons to expect that with a suitable 
redefinition of parameters the above eigenvalues T m (t\g,p)) of T m (£) are identical to the 
eigenvalues T s {t\f3 2 ,p) of the operators (|8.81| ) acting in the Virasoro module Va(c) with p 



and c already defined in (|8.78 ). The exact correspondence is 

T s (t\(3 2 ,p) = T s ( Pl t\p 2 /2,p), 8 = 0,1,2,... (8.87) 

where the factor 



Pi 



r(i-4-) r(i + l) ra-/? 2 



r(i + f ) P(i , - 

is calculated from the comparison of fl3.84|) with (|4.84j) . Further, the subspace V sym should 



support the action of the Virasoro algebra with the central charge c. It would be inter- 
esting to obtain an explicit realization of this Virasoro algebra in terms of the generators 
of the W3 algebra. It would also be interesting to obtain an explicit construction of the 
Q operators for the U q (A y 2 ; )-related T operators ( |S.81| ). 

8.3 Concluding remarks 

The study of the integrable structure of CFT in [0-f| and in this paper once again demon- 
strates the fundamental role of Baxter's Q-operators in the theory of integrable quantum 



systems. Some further applications of Q-operators in CFT can be found in ||69|| . 

Here we adopted the approach of ||[4| to construct the Q-operators as traces of certain 
quantum monodromy matrices associated with infinite-dimensional representations of the 
g-oscillator algebra. As remarked in H this construction is rather general and can be 
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applied to the lattice models as well. In particular, as the structure of the functional re- 
lations considered in Section 5 is completely determined by the decomposition properties 
of the products of representation of the quantum affine algebra U q (sl(3)), all these func- 
tional relations remain valid for the associated lattice models, with minor modifications 
related to the normalisation conventions of the lattice transfer matrices. Note that for 
certain lattice models there exists an alternative approach to Q-operators JHJ (based on 
Baxter's original work |7l]]) which allows the explicit calculation of their matrix elements 



in a simple product form (see also |72 H78|1 for further examples). It would be interesting 
to understand a relationship between the above two approaches to the Q-operators. 

The conjecture ( |7.66| ) of |2(| is a generalisation of an earlier conjecture of the same 
authors J79j which established an exact relation between the spectral determinants of the 
(second order) Schrodinger operator with homogeneous potentials |50| , |5I| and the vacuum 
eigenvalues of the Q-operators in the c < 1 CFT ||. This earlier conjecture of [79] has 
been extended for the more general Schrodinger operator (which includes the "angular 
momentum" term) 

1(1 + 1) 



d 2 x V(x) + { E - x 2a - }v(x 







(8.89) 



and proved in | 62| . Here the parameters I and a are related to the Virasoro vacuum 
parameter p and central charge c defined in fl8.78|) as I = 2p//3 2 — 1/2 and a = 1/ (3 2 — 1. 



The most essential ingredients of this proof can be summarized as follows: 
(i) The analytic properties (both in the spectral parameter E and the angular momentum 



I). For the spectral determinants of fl3.89|) these properties are deduced from the 
standard theory of the Schrodinger equations (see, e.g., |8l^|). On the other hand, 
in H and |30| the same analytic properties for the vacuum eigenvalues of the Q- 
operators were establish in c < 1 CFT. Note that the analysis of these analytic 
properties in the variable p in |30| (corresponding to I in ( |3.89| )) was based on the 



"real time perturbation theory" results (analogous to the representations ( |7.26|) of 
this paper). 

(ii) The quasiclassical approximation. The same real time perturbation theory results 
were used in [[3(J to obtain the leading large p asymptotics of the coefficient of the 
expansion of the eigenvalues of the Q-operators in the spectral parameter. These 



asymptotics exactly matched |62[ the corresponding asymptotics for the spectral 
determinants of ( |3.89|) calculated from the standard quasiclassical approximation to 
( |3.89| ) for large values of 
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(in) Quantum Wronskian relation and the Riemann-Hilbert problem. The final step of the 



proof of [62] was based on the key observation of [79] that the functional relation for 
spectral determinants of the Schrodinger equation obtained in |8(||81| (see also [|62 



for the case of general /) exactly coincide with the quantum Wronskian relation for 
the eigenvalues of the Q-operators in c < 1 CFT obtained in 0. In fact, it was 
shown |3(| that the above properties (i) and (ii) transform this functional relation 
to a simple Riemann-Hilbert problem which has a unique solution. 

It should be noted that no additional assumptions (like the phase assignment in the Bethe 
Ansatz equations for vacuum eigenvalues, similar to (p.4|)) were used in |62]. 



It would be interesting to see whether is it possible to modify the above arguments to 
prove the conjecture (|7.66|) related to the third order differential equation (|7.58|) . It seems 
that there no difficulties in generalising the above properties (i) and (ii). In particular, we 
have verified that the asymptotics of the connection coefficients Ri(E) in ( |7.65| ) obtained 
from the quasiclassical approximation to the differential equation (|7.58|) for large values of 
ii exactly match the p — > oo asymptotics of the eigenvalues of the Q-operators obtained 
from (|6.9| ) and ( |6.30|) . It is also easy to prove that the connection coefficients Ri(E)/Ri(0) 
satisfy the same quantum Wronskian relation ( |1.13 ). However, it is currently not clear 
how transform this relation to a tractable Riemann-Hilbert problem which now should 
involve two complex "angular momentum" variables. 
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Appendix A. Universal i?-matrix 



Here we give a short proof of Proposition 1, which allows us to derive the basic 
properties of the universal L-operator (|3.34|) . 

Due to ( |3.6| ) and (|3.14| ) the basic property ( |3.9| ) of the universal i?-matrix implies the 
following relations: 



(1 ® Vi + Vi ® q- hi )TZ = TZ{1 ®yi + yi® q hi ), 



(A.l) 
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that is 

[l®Vi,n\ = Kfa <g> q hi ) - ( yi <8> q- hi )K. (A.2) 
Sum up ( |A.2[ ) over i and put / = J2i Hi- We have 

3 3 

[7, 7^] = lz(J2 y l ®q K )- (J2 Vi ® ( A - 3 ) 

i=l i=l 

Let us imagine now that we check equality (|A.3|) directly. For this we compute the 
commutator on the left hand side using the relations 

[Vi, x o\ = ki q _g-i ■ ( A ' 4 ) 

After application of the Leibnitz rule we have certain series over yi ® 1 and 1 <S> %i with 
q hi and q~ hl inserted somewhere in the right tensor space. 

We move all q~ hi to the left and all q hi to the right using the relations 

q hl Xj = q " ■'•//'' . q~ hi Xj = q" * A q '' ■ (A.5) 

The result should coincide with right hand side of ( |A.3|) . 

Suppose now that we have some /, Aj and B, L , satisfying the conditions ( |3.16D , ( ^.17| ) 
of Proposition 1, and we try to compute [I (g) 1,7?.]. Again, we first use (|3.16|) , to get 
monomials over yi Cg> 1 and 1 ® x, with Ai and Bi inserted in the second tensor component, 
then move Ai to the left and Bi to the right by means of (|3.17| ). Clearly, we should get 
the same result as in the first calculation with q hi ® 1 replaced by — Bi ® 1 and q~ hi ® 1 
replaced by — Ai ® 1, since we used identical commutation relations. This proves ( |3.1S| ). 



Appendix B. Functional relations 



In this appendix we derive the determinant formulas (|5.1| ), (|5.2|) for the operators 
T M (t), and T M (t) and the functional relations (|5.4j ), (|5.5|) . The main ingredient of the 
proof is similar to that of j|] and consists of decomposition of the tensor product of 
certain degenerated Verma modules into a sum of shifted evaluation Verma modules. 

Let 7i\)~(t) be the evaluation Verma module of the algebra B + with highest weight 
ji = (fix, fi2, /U3) with respect to evaluation map ( |4.15| ). We fix the basis ?%, where n e Z> , 
that is n = (n 1; n 2 , 113) and rii > 0, rii G Z, which is given by the action of the generators 
F a ,Ff3 and F a+ p of U q (gl(3)): 

vn = t^^q-^F^FZl^vo, (B.l) 
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where vq is the vacuum vector. Then the action of the generators of the algebra B + is 
given by the relations: 



VlVn = t*[n 3 ] q \p,p - U 3 + l} q Vn-e 3 + U q^ 2n '- i [n 2 ] q Vn+e 1 -e 2 , 

1lr.1l- — /7 n 3+Ail+M3 + |,,_ 

VzVn = t^[ni\ q [^ a - U x - Tl 2 + 71 3 + l]^-^ ~ ^l?™ 2 "™ 3 " 2 "^ M^n-^+es' 

(B.2) 

hiVn = (At/3 + n 1 -n 2 - 2n 3 )t%, 
/i2fn = (~AW/3 + "i + 2n 2 + n 3 )i%, 

h 3 Vn = {fJ-a - 2ni - n 2 + n 3 )l%, 

where 

(B.3) 

ex = (1,0,0), e 2 = (0,1,0), e 3 = (0,0,1). 

Any weight \x defines a shift automorphism of the algebra B + : 

PuiVi) = Vi, Vn{hi) = ^1-^2 + ^3, , . 

(B.4) 

P/X^) = ^2 + ^1-^3, P/X^) = ^2-^1+^2- 

For a representation 7r of £> + denote by 7r[/x] the shifted representation 7r-p M . In particular, 
denote by 7f+(t) the shifted evaluation Verma module 7r+ (t) [/i] . In this module the action 
of generators is given by ( B.2|) while the vacuum vector is of zero weight and 

hiv w = (ni - n 2 - 2n 3 )vn, h 2 v w = (m + 2n 2 + 713)1%, 

(B.5) 

^3^n = (n 3 - 2m - n 2 )t%, 

so the traces are related by 

Tr nMt) e^ Fh C = tftftfTtt+Q e^C. (B.6) 

Let us define three degenerations 7f*(t), i = 1,2,3 of shifted evaluation Verma module 
TTp(t). They correspond to the limit q — > +00 in three different asymptotic zones 

Ai a > 0, /i/3 > 0; A*a ^ 0, /i/3 > and /i Q <C 0, /ia < (B.7) 
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under assumptions that all the numbers are real. 

All the degenerated Verma modules fc l (t) are equipped with bases ?%, Hi e Z> , on 
which the Cartan generators act as in ( |B.5| ) and 

V\Vn = {q-q~ 1 )~ 1 <f~ n3+1 [n3] q Vn-e 3 +g 1 '" 2n3 K]qWw+£ 1 -e 2 , 

V2Vft = q n3 + ^Vn+s 2 , ( B - 8 ) 

y 3 vn = (g-?- 1 )-V" ni " na+n3+1 [m] 9 «n-? 1 , 

for the module vr 1 (t), where we use the notation t = q 2x ; 

ViVn = (? - f 1 )" 1 ^ n3+1 [n 3 ] ? t%_ 53 + g 1 - 2 " 3 ^]^-^, 

V2Vn = q m+ ^Vn + e 2 , (B-9) 
y 3 Vn = -(Q ~ Q~ 1 )~ 1 <f+ ni+na ~ na ~ 1 [ni]qVB^i - q X+n2 - n3 - 2 [n 2 } q Vn-e 2 +e 3 , 

for the module 7f 2 (t) and 



y\Vn = ~(q-q x ) V +na H^]?^- 



£3 5 



2/2*% = q m+1 *Vn + s 2 , ( B - 10 ) 



for the module ff 3 (t) . 

The determinant formula ( |5.1| ) is based on the decomposition of the tensor product 



vf 1 ^ 1 ) <g> n 2 (tq 2 ^) ® ff 2 (V 3 ). (B.ll) 

The claim is that this tensor product admits a filtration, labelled by six arbitrary non- 
negative integers ki,k 2 , /i,/2, mi,m 2 with factormodules, isomorphic to 



it, 



'(t) [(h + k 2 )a + (h + l 2 )P + (mi + m 2 ) (a + 



where a = (1, — 1, 0), (3 = (0, 1, —1) are simple roots for the Lie algebra sl(3). 
On the level of characters this means an equality of the form 

7f 1 (tg 2 ^ 1 ) <g> Tx 2 (tq 2 ^ 2 ) ® Ti 2 (tq 2 ^) = 
T^kaum^o [(*i + fc 2 + »i + »2)« + (*i + k + m 1 + m 2 )/3], 
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(B.12) 



or, if we put Aj(£) = Tr^im e mFh £, which converges in a region ( 4.48|) , 

Mtq 2 nMtq 2 nMtq 2tls ) = A~X+(t) e^ Ph £, (B.13) 

where 

A = (1 - z 2 / Zl ){l - z 3 / Zl ){l - z 3 /z 2 ). (B.14) 

We construct the decomposition of tensor product (|B.11| ) in two steps. First, consider the 
second order tensor product vr 1 (ti) <8> vr 2 ^), where we choose the basis v™, n,m G Z> 

«?= e e (r) (?) ^™ , ( i i/^) , " /2i '«®%- E . (B.i5) 

fc 1= 0fc 2 =0 ^ X / <3 ^ 2 / ? 

where = (fci, fc 2 , 0), 

9 2 " - 1 Ai\ (n)„! 
W« = ^T> UJ, = W.I(n -,).! ■ <R16) 

and the phase Q(n,m, k) is 

Q(n,m,k) = (ni - fci)fc 2 + (&2 - fci)w 3 + £; 2 m 2 . (B.17) 

Define the following partial order in Z> : m' < m if 

m' 3 < m 3 , m' 2 + m' 3 < m 2 + m 3 and m'i + m' 2 + m' 3 < mi + m 2 + m 3 . (B.18) 

Let y m be £>+ submodules of 7r 1 (ti) ® vr 2 ^), generated by the free action of the jji, 
i = 1,2,3, 

V m = J2 (B.19) 

m'<m 

Then the vectors v^ 1 , n e Z> form a basis of the factormodule 

ym = ymj V m ' (B.20) 

m'<m 
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and the action of the generators of £>+ in V m is given by 



Vl v™ = (g-g- 1 )-V 2_n8+1 [ w 3]« u «-B8 +q X2 ~ 2n3 [n 2 } *' m 



e 3 1 H l' L 2 }q V n+e 1 -e 2 ^ 



y 2 v™ = q n -M v ™ + _ 2) 

V?,v~ = q Xl [xx - x 2 - ni - n 2 + n 3 + 

-q X2+n2 ~ n3 ~ 2 [n 2 ] q vf_ W2+ ^ (B.21) 

= ( n i _ n 2 - 2n 3 + mi — m 2 — 2m 3 )u™, 

= ("i + 2n 2 + n 3 + mi + 2?ti 2 + m 3 )v%, 
^™ = (— 2ni - n 2 + n 3 - 2m x - m 2 + m 3 )v^. 



where ti = q 2xi , t 2 = q 2x2 . 

Let 7r 12 (ii,£ 2 ) be the following representation of B + with a basis ?%, n G Z' 



>o- 



£2' 



Z/3^n = - x 2 + 1 - n! - n 2 + n3] 5 [ni],%- ?1 - <f 2 +™ 2 n ^ 2 [n 2 ] 9 u s _£ 2+ £ 3 , 

(B.22) 

(the Cartan generators act as in (|B.5|)). Then, the relations (|B.21|) mean that the tensor 



product ® n 2 {t 2 ) admits a filtration with factormodules isomorphic to 

^ 12 (h,t 2 )[(m 1 + m 2 )a + (m 2 + m 3 )/3], m, > 

and 

Ai(*i)A 2 (t 2 ) = A _1 Tr^ 12(tlit2) e™ Ph £. (B.23) 

Next, we decompose the tensor product 7f 1,2 (ti, t 2 ) ® tt 3 ^). Here we use the basis v%, 

3 

>oj 

min(7ii,n 3 ) n 2 n 3 -j / \ / \ / \ / \ 

j= o fc2 =ofc 3 =o VVgV % J q \3j q \3j q 

(B.24) 

X( ^(«) (t 3 /t 2 )^ 2 ^_^ 3) ® 
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n, m E Z 3 



where k = (0, k 2} k 3 ) and the phase £lj(n, m, k) is 

Qj(n, m, k) = j(n 2 - k 2 - n 3 + k 3 - mi - 2m 2 - m 3 ) + k 2 (k 3 — n 3 — mi — 2m 2 + m 3 ) 
+ni(2mi + m 2 - m 3 ) + n 2 {m x + 2m 2 ) + n 3 (m 2 - m x + 2m 3 ) + 3j(j 2 +1) . 

In this case we define the another partial order in Z> : m! < m if 

rn'i < mi, rn'i + m' 2 < mi + m 2 and rn'i + m! 2 + m/ 3 < mi + m 2 + m 3 . 

Let V m be the following £> + submodules of 7r 12 (ti,t 2 ) (g> vr 3 (t 3 ): 

V m =Y J B + {v ®v m ,). (B.25) 

m'<m 

The vectors v™, n e Z> form a basis of the factormodule 

ym = ymj V m ' (B.26) 

m'<m 

in which 



yi vg = q^[ X2 - X3 - n 3 + l} q [n 3 } q v^ 3 + q X2 2n3 [n 2 ] q v^ +Ii _ l2 , 



y 2V ™ = q ^+k v ™ + ^ 



y 3 v% = q Xl [x 1 - x 2 - m - n 2 + n 3 + l] g M g ^_ £l - 



(B.27) 



_„X2+n 2 -n 3 -2r 1 m 



with tj = g 2Xl and the Cartan generators acting as in (p.21| ). We see that the factormod- 
ules are shifts of the same £> + -module which is isomorphic to vr^(t), if we put £, = tq 2lli . 
This in combination with ( |B.23| ), proves the main statements ( |B.12| ) and ( |B.13| ). 



In order to deduce ( |5.1|) from (|B.13|) we notice that the degenerated evaluation Verma 



modules 7r l (tj) are not irreducible; each of them admits a filtration indexed by a single 
nonnegative integer such that all the factors are the shifts of oscillator representations: 
the representation 7f 1 (t) admits an increasing filtration by submodules V n , where V n 
is generated over B + by the vectors i>(o,o,fc)> k < n with factormodules isomorphic to 
representation pi(q 2 t)[n/3] ,n > 0; the representation 7r 3 (t) admits an increasing filtration 
by submodules V n , which are generated over B + by the vectors vtk,o,o), k < n with 
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factormodules isomorphic to representation p 3 (q~ 2 t)[na] ,n > and the representation 
Tf 2 {t) admits a decreasing filtration by submodules V n , which are generated over B + by 
the vector 

gwi^^fe-'-^W.. (B.28, 

The factormodules are isomorphic to p 3 (t)[n(a + /3)], n > 0. 
This means that 

Ax(t) = (l-z 3 /z 2 )- 1 Tr pi{q 2 t) e™ Ph C = A~ 1 A 1 (qH), 
A 2 (t) = (l-z 3 M)- 1 Tr p2(t) e- Ph £ = A^A 2 (t), (B.29) 
A 3 (t) = (1- z 2 /z 1 r 1 Tr ps{q - H) e i ^ h C = A" 1 A 3 (q- 2 t) . 
The relation (|B.13| ) combined with ( B.29 ) gives 

Tr-+ (t) e lwPh C = A~ 1 A 1 (tq 2fll+2 )A 2 (tq 2 < l2 )A 3 (tq 2 ^- 2 ). (B.30) 

The determinant formula (|5.1| ) is a direct corollary of ( |B.30| ) and of the Bernstein- Gel'fand- 
Gel'fand resolution of finite-dimensional module n^, which implies the relation 

Tm(«) = E(- 1 ) ,Wt 2Wp)-p(*)> ( R31 ) 

<t£S 3 

where p = (1, 0, —1). 

The determinant formula (|5.2|) can be deduced in an analogous manner with the use 
of the second evaluation map or can be deduced from ( |5.1| ) in the following way. Due to 
the definition (|4.18j) , we conjugate the action of the generators of B + in the shifted Verma 
module nt(t) by cr 13 . Taking into account the fact that the U q (gl(3))-Verm& module 
ft 111,112,1*3 ' °"i3 * s isomorphic to the Verma module it-^-^ _ w , we see that 

%(t) = ^ i3wH . (lte) e™ Fh C and T+(f) = (-*).(«*!.) e^ Ph £, 

which means that we can use the statement of ( |5.1|) for the highest weight — <x 13 (/x) com- 
bined with automorphism a\ 3 of the algebra £>+. Due to the definition ( |4.46| ) we get the 
statement of ( |5.2j) . 

Let us turn to the relations ( |5.4|) , (|5.5| ). Their proof consists of an interpretation of 
both sides of the equalities as of £> + -modules. Consider, for instance, the tensor product 

n 2 (t 2 )®n 3 (t 3 ), 
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of degenerated evaluation Verma modules. Analogously to ( |B.22|) , ( |B.23|) it is isomorphic, 
on the level of characters, to a sum 

22 K 23 {t 2 , h) [(mi + m 2 )a + (m 2 + m 3 )P)), 

m,i>0 

where the module vr 23 (t2, H) is the following representation of £>+ with a basis n G Z> : 



ViVn = <f z [x 2 -x 3 -n 3 + l} q \n 3 } q Vn-e 3 + q X2 ~ 2n3 [n 2 } q Vn + - 
V3Vn = -(q - q~ 1 )~*q"*~ r " ,1 ~ r ""-"'-*[n 1 \ q Vn-- Sl - </ 



£1—62 1 



(B.32) 



where hi act as in ( |B.5|) , t; L = q 2xi . When t 2 = t 3 = t, n 23 (t 2 , t 3 ) has a submodule generated 
by the vector Ve 3 with a basis ■% = fn+e 3 , n G Z> . One can see from ( |B.33| ), ( |B.5|) that 
this submodule is isomorphic to 7r 23 (q~ 2 t,q 2 t)[/3]. The factormodule 

^(t,t)/n 23 (q- 2 t,q 2 t)[P} 

has a basis f n m , n, m > with the action of £> + given by 



yiv n ,m = q x ~ 2p [m)qV n +i,m-i, 

y2v n ,m = q*v njm+1 , 

y?,v n ,m = ~(q - q~ 1 )~ 1 <f~ H *™~ 1 [n] q v n - ltm , 

Mn,m = (n- m)v n ^ h 2 v n , m = (n + 2m)v n , m h 3 v n>m = (-2n - m)v ri 



(B.33) 



From ( B.33| ) we can check directly that this factormodule is isomorphic to the represen- 
tation pi (q~H), so 

A ( A 2 (t)A 3 (t) - z 3 /z 2 A 2 (q- 2 t)A 3 (q 2 t)) = (1 - z 3 / Zl )-\l - z 2 /z 1 )- 1 A 1 (q-H), 



or 



A 2 (gt)A 3 (g- 1 t) - z 3 /z 2 A 2 (q-H)A 3 (qt) = (1 - z 3 /z 2 )A 1 (t), 



(B.34) 



which is equivalent to one of the relations (|5.4j) . The others follow from (P-34Q after 
application of the symmetries (|4.9|) and (|4.10| ). 
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Appendix C. The p — > oo asymptotics of ^(V). 

In this Appendix we derive the leading asymptotics of the function *f?(v) at large 
values of p = (pi,p2)- We assume pi and P2 tend to infinity while their ratio P2/P1 is kept 
fixed 

p 1 ,p 2 -> 00, p 2 /Pi = finite. (C.l) 
First consider the classical case. Writing the formula (|2.28 ) in full one gets 



T^(A) = Ai d) (A) + Af \X) + A 3 C °(A), Kf\X)A.f\x)A.f\X) = 1, (C.2) 

where Aj c ^(A), j = 1,2,3 are the eigenvalues of the classical monodromy matrix M(A). 
The leading p — > 00 contribution to these eigenvalues can be easily calculated directly 
from the definition fl2.27| ). According to ( |2.44| ) the potentials U(u) and V(u) in ( |2.26| ) in 
this limit become independent of u 

^(M)|p^oo ~ XiX 2 + X1X3 + 2:2X3, V(m)|p^oo ~ ixix 2 x 3 , (C.3) 

where 



x 



j 



[w 3 p), j = 1,2,3 (C.4) 



with m,- defined in (|2.34[ ). Therefore the eigenvalues of M(A) 

Kf\\) = e 2mx ^ x \ (C.5) 
are expressed though the roots Xj(X), j = 1, 2, 3, of the cubic equation 

x(X) 3 + g 2 x(X)+g 3 -iX 3 = Q, (C.6) 
where we have introduced the notation 

g 2 = x x x 2 + X1X3 + x 2 x 3 , g 3 = -xix 2 xi. (C.7) 
The roots of (|C.6|) can be written as 

m 9IY 92\h ( 111, 27^3-^A 3 ) 2 ,] 

s i A =2 (- — ) 2 ^i(--, -,-;! + j-2 ) > ( c - 8 ) 

LV 3 / o 2 4^ -I i- th 

where 2 Fi stands for the Gauss hypergeometric function. The expression in the square 
brackets above is a multivalued function and the symbol denotes an appropriate 

branch of this function fixed by the requirement Xj(0) = Xj where Xj, j = 1,2,3 defined 
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Ill 



Expanding this expression in a series in A and using the Sommerfeld- Watson 

as 



transformation one can rewrite QC.8 



logAj 
where 



iv \ 



(A) = 2^ + I -r(i - ^)r(-^ + "-)rQ + ^-)(-zA)t *f (*,), (c . 9 ) 



7/ 



2"^ 



37iT(l/2 + ii//3) 
# 2 \ (l-i^/2 



x 



3 / 



1 ii/ 1 iz/ 1 iz/ 27(7 
l( ~6 + 6"' 6 + 6"' 2 + 3"' 1 + 1^1 



J 7-th 



rg.io) 



where the integration contour goes along the line 3m v = — 1 — e and the branches of the 
expression in the square brackets are chosen in accordance with ( |C8| ) 



*; d) (o) 



WjP) 



(C.11) 



27^' 

Now turn to the quantum case. As before we restrict our considerations to the case 



< g < |. Examining expressions (|6.9|) it is easy to conclude that the leading p — ► oo 
asymptotics of ^(f) is the same for all its eigenvalues, i.e. ^(z^lp^oo is proportional to 
the unit operator I. Further we make an assumption (which again can be justified by 
studying the p — > oo of the Coulomb integrals (|4.75|) ) that the p and g dependence of 



9(u) 



factorizes in the form 



Comparing ( |C.9| ) with the g — > limit of ( p.9| ) one gets 



^ 9) -*f\v)\ P - 



(C.12) 



(C.13) 



To calculate the (^-dependent multiplier in the last formula we will use the DDV equation. 
Since the asymptotics (|6.16|) are the same for all eigenvalues of *&(v) it is enough to 
consider the vacuum eigenvalue. Moreover, since the multiplier we want to calculate does 
not depend on p we will restrict ourselves to the case P2 = and denote 



Pi =P, 



x 1 = p, 



x 2 = 0, 



x 3 

r(vac) 



-p- 



(C.14) 



In this case the vacuum eigenvalues A^^it) and A£"^'(t) of the operator Aj(t) and 
As(t) satisfy the relation 

(C.15) 



.4 



(vac) 



(t) = A^ ac \-t). 
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Define the function 

A ( r C \itq- 2 ) A { ^\-itq) 
Recalling the assumption about the structure of the zeros of A\(i) in the beginning of 
Section ^ (the zeroes of (i) lie on the positive imaginary axis of t) and using ( |6.7| ) 
we can show for real p 

a (t)* = a^*)- 1 . (C.17) 
Due to ( |6.1| ) a{t) has the asymptotics 

log a(t) 2zMsin ( j^^^ J t ^ oo, | arg(t) | < min[27r#, 7r(l - g)} 

(C.18) 

while for small t, we obviously have 

a(f) = 2mp + 0{t), t^O. (C.19) 
Again using ( |6.7|) we have 

log a(t) = 2mp + F ^ u k l ), (C- 20 ) 

k 

where 



n«) = ^ ( t^^- ) ■ (c.2i) 



1 - tq 2 1 +tq- r 
l-tq~ 2 T+tq 

Using the Bethe Ansatz equations (|6.3| ), (|6.4j) we can write the sum in ( C.20 ) as a contour 
integral 

\oga(t) = 2mp + / -^Fitv-^dylogil + aiv)), (C.22) 
Jc 2m 

where the contour C goes from +oo to zero above the real axis, winds about zero, and 
returns to +oo below the real axis. Integrating by parts (boundary terms vanish) one 
obtains 

-j POO I 

log a{t) - 2mp =— — tdtFitv- 1 ) {log(l + a(v - iO)) - log(l + a(v + iO))} , 
2m Jo v 

r°° dv f°° dv 

= —td t F(tv- 1 )$hnlog(l + a(v-iQ))- td t F(tv~ l ) loga(u), 

Jo J 2mv 

where we have used ( |C.17| ) to obtain the second line. 

We now introduce new variables 

39 39' 

t = e~i, v = e~i, (C.23) 
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and recalling that q = exp(ing), we obtain 

/oo 
R{6 - 6')^m\og{l + a{6' - zO)) 6ff 
■oo 

/oo 
R{B-&)\o%a{&)dff, (C.24) 
-oo 

where 

#(0) = — - — -td t F(t). (C.25) 

With the standard notation for the convolution integral the Eq. QC.24 ) becomes 

K*\oga(9) =2irip-2iR*5imlog(l + a(6 -iO)), (C.26) 

with i^(0) = 5(8) — R(6). We now apply the inverse of the integral operator K(9) to 
both sides of (|C.26|) . Since K{6) is a singular operator one has to add an appropriate zero 
mode of K to the RHS of the resulting equation. The required zero mode is completely 
fixed by the asymptotic conditions ( |C.18| ) and (|C.19|) . In this way one obtains 

i log a(6) = - 2G * 3m log(l + a(8 - iO)) + 2M sin y ^^ J e °> ( C ' 27 ) 

where G(0) = -(K' 1 * R)(6). 
Now define 

B(p) = 2-^tog(-ft 1 ), (C.28) 

where t\ is the smallest (in absolute value) zero of Al ac {t). The convolution term in (|02l ) 
can be written as 

fB( P ) 

G*3mlog(l + a wc (0-iO)) = / G(0 - 6')Smlog{l + a vac {6' - iO)) dB' 

J —oo 

POO 

+ / G(6 - 6')%m\og(l + a vac (6' -i0))d6'. (C.29) 
Jb(p) 

Let us make the technical assumption that 

B(p) ~ const log p, p — > oo. (C.30) 

Then one can show the second term in (|C.29| ) vanishes when p — > +00 and therefore can 
be dropped in the leading large p approximation. The Eq. ( |C.27|) then becomes 

_ HE + Msin f^L^I e e = / B(P) ^d \og S(0 - 003m log (1 + a(6')) , (C.31) 
9 V 3 / ./-oo 2m 
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where we have used the notation 



G(6) = 6(6) + ^d e log S(8), (C.32) 



sinh ( cosh ) 
and the relation 

\oga vac (6) = 2zlog3mlog(l + a mc (#)), 6 < B{p) (C.34) 

which holds for 6 < B(p). The function S(6) is Smirnov's kink-kink S'-matrix for the 
algebra |)6 . 



The Fourier transform of the integral kernel in (|C.31|) factorizes as 

i r°° i 

- 7T- / e-™ e de log S{9)dB = (C.35) 



where 



i r (-s^fl) r (| - f ) e< 



= ^ r(- ¥ )r(i- f ) ■ *-<") = *+<-'>■ < c ' 36 » 

3<S = (1+0 log(l + - { log { + log (C.37) 

The function K+(y) is analytic and nonzero in the upper half-plane 3m u > — 1 where it 
has the asymptotics 

K+(v) ~ 1 + O ( - J , zv^oo. (C.38) 



The Eq. (|C.31|) is an equation of Wiener- Hopf type, which can be solved by the standard 
technique. Defining the Fourier transform 

fB( P ) 

f+{v)= ^m\og{l + a vac {6))e- iev d9, (C.39) 

J — oo 

one obtains 

lp^L e iuB ^ = — -^— K (v')G + (u')du', (C.40) 



where 



inp %M { ir{l + £)\ _ B(p) 



gz/ 1/ + i y 3 
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sin -^-^ ^W, (C.41) 



and e is a small positive constant. Completing the contour in the lower half-plane, we 
have 

AM>%) = ™2_ _ iM™ B <*)- s r (|) 

tf+M " ^ ^ + r(f)r(^)' 1 ' ] 

As follows from the definition ( |C.39|) the leading \jv term in the v — > oo asymptotic 
expansion of the RHS of (|C.42| ) is independent of p. This gives 

(-^i)^= / ( |r ( y ) e^ (C.43) 



and thus we have 



^p-^jaV ^^- (a44) 



Given the function f+(v), QU.44 ), we can recover log Ai (t) from 



where the contour Cj, is along the line Sm 1/ = — 1 — e, with e is a small positive constant. 
Thus we have 



^Ar{t) = -^-l e —^ — ^W-^ — -(—) 



7c, ^ 2 r (| + f ) 

(C.46) 

This expression can be written in the form of (|6.9| ) 



i()gA,(-/-) = — / — r(K/3)r(i-iKe + i)/3)r((-i + ^)/3) x 



27ri 7(7^ 



r((i + .v)/3) (r(i - (it)^ 1 -^/ 3 ^WUoe, (c.47) 

where 



2-TT / p 2 N (l-^)/2 



* (iv)l ^°°" V3 l^J r(i/3 + ii//6)r(2/3 + ii//6)- (C-1 '' 



Comparing this with the specialization of (|C.10|) for the case flC.14p , one concludes that 
the factor in ( |C.13| ) 

#4 es 1. (C.49) 
Taking into account ( |C.13| ), we obtain the expression ( |6.30| ) given in the main text. 
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To get a small t expansion of A^ ac (t), we complete the contour in the lower-half plane, 
and evaluate the residues at the points v = 3i(g — l)n, n = 1, 2, 3 . . ., and we obtain the 
large p asymptotic equivalent of ( [4.55 ). Thus, we have from flC.46| ) 

l P ^oo- 2r3(1 _^ r(i-i^) M l} • ( } 

However, we know /f| 1 ^' uoc ' exactly ( [4.851) . Taking the large p asymptotic of (|4.85| ), 
comparing with ( |C.50| ), and using ( |C.43j ) we see that we must take 

M = M £ 3 > (r(l - g))~ , (C.51) 

and then 

r-iv*r( 3(1 -f- l N )r(gn) 

#«( vac ) L^oo = J — s ^3y- 3(1 - 9)n - (C52) 

lp 2n\ r Ai-g)"+i x 



2 

We also generate an asymptotic expansion for large t by completing the contour in 
the upper half-plane. Here we have poles at the following points: 

i/ = i(2n-l), n = 0,1,2,..., v = 3i(g- x - l)n, n = 0, 1, 2, 3 . . . , (C.53) 

and we obtain 



Ar{t) = c va \g,p)(s)-U W {Y J B 2n -JZ-i\^oo{s)^ 



3(1-9) 



. n=0 



x exp -j^fl-W^l^oo^-f , (C.54) 



n=l 



where 



-lr 1 r (*Sy) r G 



2n-l 



^- 1= 3(l-,)(n)! * 9 ' (C ' 55) 

C^\g,p) = ( H!£ ) 9 (^)-f , HV>W(g,p) = H^ ac \g~\ g- l p). (C.56) 

and I^iilp-Kxjj Ji = 0, 1, 2, . . . are the remaining Integrals of Motion after the A 2 spe- 
cialization has been taken into account (W(u) = 0), 



p2n 

jwc _ -^2n-llp^oo = - — , w = 1, 2, 3, . . . (C.57) 
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Note also that C {vac \g,p) is the vacuum asymptotic value of the operator Cj{g, g<p{u)} 
defined in (|6.16|) . 



Appendix D. Representation Matrices for the Funda- 
mental Representations 



For the first fundamental representation 7r Wl we have 



\ 



7T Wl {Hp 



1 
0-10 
\ ) 

( o o o N 

1 
0-1 



/ 



^{Eo 



while for 7r W2 we have 



n W2 {H Q ) 



^ {Hp 





1 

0-1 

1 
0-10 




X 



wi {Ep 



7^2 {E a 



«u z {Ep) 
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y j 

( o o o N 

1 




x^{Ea, 
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1 



1 
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X^) 2 { Ea j 



K0J2 {Fp 
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(d.i; 



(D.2) 



(D.3) 



(D.4) 
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